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1.0  Final  Report  Summary 
Technical  Problem 


The  technical  problem  addressed  by  this  contract  was  to 
determine  the  feasibility  of  exploiting  the  emerging  mathe¬ 
matical  theory  of  wavelets  in  radar  system  applications. 

The  work  to  be  done  was  divided  into  three  tasks. 

1.  Perform  a  preliminary  systems  benefit  study.  Address 
the  complete  radar  processing  problem  including  the 
applicability  of  wavelets  to  the  various  processing 
elements:  prescreener,  fast  Fourier  transform  (FFT), 
classifier,  and  automatic  target  recognizer  (ATR). 

2.  Define  and  develop  appropriate  wavelet  methods  for  radar 
appl i cations. 

3.  Determine  the  stability  of  the  wavelet  transform  image. 
Make  a  comparison  of  wavelet  transforms  and  FFT  signatures 
for  each  of  three  targets:  a  set  of  calibrated  corner 
reflectors  of  known  radar  cross  section  and  spacing,  an 
M60  tank,  and  an  M35  truck. 

Three  specific  potential  applications  of  wavelets  methods  to 
radar  systems  were  identified  for  investigation. 

1 .  Use  wavelets  to  transform  the  data  that  forms  the  input  to 
the  radar’s  ATR  algorithms. 

2.  Modify  a  typical  radar  system  so  that  if  a  wavelet 
transform  (which  is  faster  than  an  FFT)  replaces  the  FFT 
(which  is  currently  used),  the  output  is  still  a  high 
range  resolution  target  profile. 

3.  Develop  a  wavelet  approximation  to  an  FFT  such  that  the 
wavelet  approximation  is  sufficiently  accurate  and  faster 
than  an  FFT. 

Methodology 

In  this  study,  the  recent  mathematical  theory  of 
wavelets  was  introduced  to  the  engineering  problems  of 
designing  radar  systems,  radar  processors,  and  radar 
algorithms.  The  goal  was  to  make  radars  more  efficient  or 
more  effective  by  the  use  of  wavelets.  To  understand  why 
particular  possible  applications  of  wavelets  to  radars  were 
examined,  it  is  necessary  to  understand  some  background 
information  on  both  radars  and  wavelets  theory.  These  topics 
are  discussed  briefly  in  the  following  paragraphs.  Also, 
the  Martin  Marietta  radar  data  that  was  used  in  this  study 
is  described,  and  some  description  is  given  of  typical 
processing  methods  for  this  data. 

Modern  radar  systems  transmit  and  receive  frequency- 
stepped  radio  frequency  (RF)  waveforms.  The  radar  signal 
processor  transforms  the  received  data  by  using  an  FFT. 

Since  the  FFT  forms  a  matched  filter  with  the  frequency- 
stepped  waveforms,  the  vector  whose  components  are  the 
magnitudes  of  the  complex  data  output  from  the  FFT  forms  a 
high  range  resolution  profile  (HRRP)  of  the  objects  in  the 
radar  beam.  These  HRRPs  are  the  inputs  to  the  radar’s  ATR 
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processor.  Both  FFT  and  ATR  calculations  are  computation 
intensive  and  impose  difficult  and  costly  requirements  on  the 
radar  signal  processor.  Wavelet  methods  were  examined  for 
potential  efficiency  improvements  in  FFT  or  ATR  processors  as 
well  as  for  potential  ATR  performance  improvements. 

Investigations  utilized  Martin  Marietta’s  short  pulse, 
frequency-stepped,  fully  calibrated  35  Ghz  radar  data.  This 
data  utilized  two  transmit  (right  hand  circular  and  left  hand 
circular)  and  two  receive  (right  hand  circular  and  left  hand 
circular)  radar  polarizations  (right  and  left  circular)  and 
resulted  in  a  complete  polarization  set  of  four  receive- 
transmit  pairs.  Transmitting  right  or  left  and  receiving  the 
opposite  yields  “odd"  polarized  returns.  If  the  receiver  has 
the  same  polarization  sense  as  the  transmitter,  the  return  is 
an  "even"  polarization.  The  data  is  coherent  which  allows 
creation  of  HRRPs.  Figures  3.3  through  3.6  show  odd  and  even 
HRRPs.  There  are  radar  looks  at  an  M60  tank  and  an  M35  truck 
at  every  hundredth  of  a  degree  of  aspect  in  this  data  set.  A 
radar  look  consists  of  64  coherent,  frequency  stepped  time 
samples  for  the  four  recei ve-transmi t  polarization  pairs. 

Each  polarization  of  each  radar  look  can  be  processed  through 
an  FFT  to  produce  an  HRRP  of  the  vehicle  in  the  radar’s  view. 
Compactly  supported  wavelets  are  a  family  of  recently 
discovered  (1986)  mathematical  functions  that  have  been 
successfully  applied  to  problems  in  image  compression,  audio 
compression,  vision  analysis,  and  transient  signal  analysis. 
The  wavelet  methods  used  on  this  contract  employed  compactly 
supported  wavelets  that  generated  orthonormal  bases.  Par¬ 
ticular  wavelets  used  were  denoted  D2,  D4,  D6 ,  and  D8  in 
honor  of  Ingrid  Daubechies,  who  developed  the  theory  of 
orthonormal  bases  of  compactly  supported  wavelets.  The  D2 
wavelet  is  also  known  as  the  Haar  wavelet.  A  wavelet  basis 
consists  of  a  "scaling  function"  0(x)  that  is  the  solution  to 
a  recursion  equation  derived  from  orthogonal  subspace 
projection  mathematics,  a  basic  wavelet  function  woo(x)  that 
is  derived  from  the  scaling  function,  and  a  countably 
infinite  collection  of  wavelet  functions  Wj k ( x )  that  are 
dilations  (frequency  changes)  and  translations  (time  shifts) 
of  the  basic  wavelet  function.  The  recursion  equation  is 
given  by 

0(x)  =  Zck0(2x-k) 

and  the  wavelet  functions  are  defined  by 
woo(x)  =  Z(-1  li'ci -k0(2x-k  ) 

Wjk(x)  =  2-1  /  2  (  2j  x-k  ) . 

Some  specific  linear  operations  involving  the  recursion 
coefficients  Ck  can  compute  wavelet  coefficients  and  wavelet 
approximations  to  functions  very  rapidly  since  only  a  few  of 
the  Ck  are  nonzero. 

Not  all  wavelets  have  compact  support  nor  do  they  all 
generate  orthogonal  bases.  However,  throughout  this  report 
the  terms  wavelets  and  compactly  supported  wavelets  are  used 
interchangeably,  and  wavelet  bases  are  assumed  to  be  orthog¬ 
onal.  Results  presented  were  produced  using  specific 
compactly  supported  wavelets. 
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Technical  Results 


Results  clearly  demonstrate  that,  under  certain  circum¬ 
stances,  wavelets  methods  can  improve  radar  system  perform¬ 
ance.  These  wavelet  results  and  their  radar  system  appli¬ 
cations  are  summarized  below  and  described  in  detail  in 
Sections  3.2,  3.3,  3.4,  and  3.5  of  this  report. 

The  radars  considered  in  this  study  were  millimeter  wave 
(MMW)  radars,  but  results  obtained  are  also  applicable  to 
radars  of  other  wave  lengths.  MMW  radars  are  typically 
either  pulse  radars  with  short  radar  pulse  lengths  or 
frequency  modulated  continuous  waveform  ( FMCW )  with  long 
radar  pulse  lengths.  There  is  a  major  difference  between 
these  two  types  of  MMW  radars  as  far  as  wavelets  applications 
are  concerned.  This  is  due  to  the  fact  that  if  one  wants  a 
given  range  resolution  in  the  HRRP,  one  must  compute  an  FFT 
of  very  large  dimension  for  an  FMCW  radar  but  a  relatively 
small  FFT  for  a  pulse  radar. 

The  most  significant  application  of  wavelet  methods  to 
radars  involves  the  radar’s  ATR  algorithms’  capability  and 
complexity.  The  data  base  of  full  polarization,  calibrated, 
coherent  MMW  radar  returns  from  tanks  and  trucks  was  used  for 
this  study.  This  data  base  also  contained  radar  looks  at  the 
calibration  reflectors  consisting  of  a  single  trihedral  and 
a  pair  of  dihedrals  of  known  radar  cross  section.  The  major 
goal  was  to  reduce  the  dimension  of  the  HRRP  radar  data 
vector  without  losing  information  that  is  important  to  the 
ATR.  There  are  three  distinct  reasons  why  this  dimension 
reduction  could  improve  the  ATR. 

1.  The  reduced  dimension  that  serves  as  the  input  to  the  ATR 
may  allow  use  of  a  more  complex  and  effective  (higher 
order  polynomial  classifier,  for  example)  classifier  or 
target-clutter  discriminator  algorithm  within  the 
computational  limitations  of  the  ATR  digital  signal 
processor . 

2.  Reduced  dimension  may  make  the  choice  of  the  type  of  ATR 
classifier  and/or  target-clutter  discriminator  obvious; 
consequently,  ATR  performance  could  be  improved. 

3.  Even  if  reduced  dimension  does  not  suggest  a  change  in 
the  type  of  classifier  or  discriminator  algorithm,  it  may 
greatly  reduce  the  number  of  computations  required 

by  the  ATR  thereby  increasing  the  speed  and  decreasing  the 
hardware  requirements  of  the  ATR. 

As  a  test  case  for  the  ability  of  wavelets  to  reduce 
HRRP  dimension,  one  feature  -  "target  length"  -  was  examined. 
For  this  effort,  Martin  Marietta  defined  the  problem  and 
monitored  the  results  while  our  subcontractor.  Aware,  Inc., 
provided  the  wavelets  expertise  and  produced  the  results. 
Wavelet  algorithms  and  typical  target  length  estimation 
techniques  were  applied  to  the  radar  HRRP  data  base.  Results 
showed  that  when  appropriate  wavelet  methods  were  applied 
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with  the  D4  or  D6  wavelet  basis,  vector  dimension  could  be 
reduced  by  a  factor  of  four  with  essentially  no  loss  in 
ability  to  estimate  target  length.  The  main  results  are 
given  in  Section  3.3  and  additional  details  are  presented  in 
Section  4.0  and  Appendix  E. 

Next  a  brief,  two  class  target  classification  study  was 
conducted  by  Martin  Marietta  to  see  if  it  is  feasible  to 
improve  radar  target  classification  signal  processing  by 
wavelet  methods,  either  by  improving  classifier  percent 
correct  classification  (PCC)  in  the  presence  of  noise  (added 
to  the  very  clean  radar  target  data)  or  by  decreasing  the 
number  of  computations  required  by  the  classifier  processor, 
or  both.  In  this  study,  the  features  used  by  the  classifier 
were  particular  wavelet  coefficients  produced  by  the  target 
HRRP.  Classifier  overall  PCC  was  improved  in  realistic 
SI gnal/interference  situations  by  up  to  6  percentage  points, 
and  computations  were  decreased  by  a  factor  of  10.  These 
results  clearly  indicate  that  it  is  feasible  to  improve  radar 
classifier  processors  by  wavelet  methods.  Details  of  this 
study  are  given  in  Section  3.4. 

A  second  application  of  wavelet  methods  to  a  radar 
system  consists  of  replacing  the  radar’s  FFT  by  a  wavelet 
transform  (WT).  This  results  in  a  more  efficient  processor 
for  two  reasons.  An  FFT  requires  on  the  order  of  NlogaCN) 
multiplies,  but  a  WT  requires  only  on  the  order  of  N 
multiplies.  Also,  an  FFT  uses  compl ex*compl ex  multiplies, 
which  requires  4  real  multiplies,  1  add,  and  1  subtract,  but 
the  waveform  proposed  for  the  WT  radar  application  is 
processed  by  a  WT  with  only  real*real  multiplies.  The  output 
of  the  WT  will  still  be  an  HRRP  if  the  radar  transmitter  and 
receiver  can  be  modified  so  that  the  received,  sampled 
waveform  data  vector  is  a  matched  filter  with  the  WT ,  just  as 
the  current  frequency  stepped  radar  data  vector  is  a  matched 
filter  with  an  FFT.  This  wavelet  radar  application  would 
improve  either  pulse  or  FMCW  radar  digital  processor  effi¬ 
ciency.  Details  of  this  WT  radar  application  and  an  example 
are  given  in  Section  3.5. 

Another  wavelet  radar  application  consists  of  using  a 
wavelet  fast  Fourier  transform  approximation  (WFFT),  derived 
by  Aware,  Inc.,  instead  of  using  an  FFT.  This  procedure  does 
not  require  any  modification  to  the  radar  transmitter  or 
receiver  but  does  require  changes  in  the  radar  digital 
processor.  A  number  of  versions  of  WFFT’s  were  examined. 

The  fastest  version  of  a  WFFT  is  faster  than  an  FFT  in 
creating  a  target  range  profile  if  the  target  length  (in 
post-FFT  cells)  is  less  than  one  fourth  of  the  FFT  dimension. 
In  cases  when  the  target  would  occupy  more  than  one  fourth  of 
the  post-FFT  cells,  the  FFT  will  be  at  least  as  fast  as  the 
WFFT.  Consequently,  the  WFFT  does  not  appear  to  be  useful 
for  typical  MMW  pulse  radars  but  would  be  useful  in  typical 
FMCW  radars  which  sample  relatively  long  waveforms  to  produce 
the  data  vectors  that  are  input  to  large  FFTs  (or  WFFTs).  A 
discussion  of  this  WFFT  vs.  FFT  comparison  follows  in 
Section  3.2.  Additional  details  are  given  in  Appendices  B.2 
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and  B.3. 


Further  Research 

This  wavelets  radar  applications  feasibility  study  has 
pointed  out  several  areas  for  further  useful  research. 

Wavelet  dimension  reduction  capabilities  may  improve 
other  radar  ATR  algorithms,  either  by  improving  stability 
versus  signal-to-interference  and/or  by  decreasing  compu¬ 
tational  complexity.  Also,  wavelets  methods  may  prove  even 
more  useful  in  two-dimensional  ATR  applications  such  as 
synthetic  array  radar,  infrared  imaging,  or  electro-optical 
i mag i ng . 

If  there  is  an  FMCW  system  whose  targets  are  small  in 
length  relative  to  the  radar  pulse  length,  this  study 
indicates  that  it  would  be  appropriate  for  that  system  to 
further  investigate  implementing  WFFTs  rather  than  FFTs. 

Several  questions  remain  concerning  the  possible  use  of 
WTs  to  replace  FFTs  in  a  radar.  What  wavelet  basis  produces 
the  radar  waveform  that  is  most  appropriate  for  transmission 
and  reception  by  a  radar  system  such  that  the  WT  would  output 
HRRPs?  Can  such  a  system  be  built?  Is  the  increased  effi¬ 
ciency  of  the  WT  (over  the  FFT)  sufficient  to  offset  the  cost 
of  the  changes  required  in  the  transmitter  and  receiver? 
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2 . 0  Compactly  Supported  Wavelets 

A  major  reason  for  conducting  this  study  is  the  remarkable 
results  being  obtained  through  the  application  of  compactly 
supported  (zero  outside  a  finite  interval)  wavelets  to  signal 
processing.  Compactly  supported  wavelets  are  a  class  of 
mathematical  functions  that  were  discovered  in  1986.  Some  of  the 
particular  advantages  of  wavelet  signal  processing  methods  are: 

o  Wavelet  transforms  are  computationally  efficient.  The 
number  of  arithmetic  operations  required  to  perform  a 
wavelet  transform  is  linearly  proportional  to  the  number  of 
input  data  points.  The  computational  complexity  of  the 
more  traditional  Fast  Fourier  Transform  (FFT)  is 
proportional  to  the  number  of  input  data  points  times  the 
logarithm  (base  2)  of  the  number  of  input  data  points.  For 
large  problems,  the  wavelet  methods  require  only  a  fraction 
of  the  number  of  operations  required  by  the  traditional 
methods.  This  advantage  increases  as  the  problem  size 
increases.  In  addition,  wavelet  transform  algorithms  can 
be  directly  implemented  in  very  large  scale  integration 
(VLSI)  logic  devices,  and  they  are  fully  parallelizable . 

o  Wavelet  transform  methods  can  analyze  signals  in  both  the 
time  and  frequency  domains.  The  relative  resolution  of  the 
time  and  frequency  components  can  be  flexibly  adapted  to 
the  problem  at  hand.  The  selection  of  the  appropriate 
time-frequency  resolution  can  be  done  upfront  at  system 
design  time  or  it  can  be  accomplished  with  real-time 
adaptive  algorithms.  The  traditional  Fourier  transform 
suffers  from  very  poor  (or  nonexistent)  time  resolution. 
This  particularly  limits  its  usefulness  in  the  analysis  of 
time-limited  (i.e.,  transient)  signals.  There  have  been 
attempts  to  modify  the  Fourier  technique  in  various  ways  to 
overcome  this  limitation,  but  all  of  the  methods  introduce 
some  additional  complexities  and  compromises.  Wavelet 
methods  offer  a  very  natural  means  to  perform 
time-frequency  signal  analysis. 

o  Wavelets  provide  the  flexibility  to  choose  a  particular 
wavelet  function  that  is  "customized"  to  the  specific 
application.  This  is  possible  since  compactly  supported 
wavelets  are  an  infinite  family  of  complete  orthogonal 
basis  functions.  This  flexibility  to  choose  basis 
functions  can  not  be  matched  with  the  Fourier  transform  for 
it  uses  only  a  single  set  of  basis  functions  -  the  complex 
exponentials  (i.e.,  the  sine  and  cosine  functions.) 

Compactly  supported  wavelets  are  a  complete  and  orthogonal 
set  of  basis  functions  for  the  set  of  all  finite  energy  discrete 
signals.  The  wavelet  transform  is  invertible,  energy-preserving 
and  linear.  The  wavelet  transform  is  a  processing  method  which 
analyzes  both  continuous  streams  of  input  data  and  blocks  of  data. 
A  multiplier  2  wavelet  basis  consists  of  a  scaling  function,  a 
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basic  wavelet  and  a  collection  of  smaller  wavelets.  The  smaller 
wavelets  are  created  by  "shrinking”  the  basic  wavelet  by  a  factors 
of  2  and  shifting  (or  translating)  them  by  scaled  integer 
distances.  Thus  the  collection  of  smaller  wavelets  are  1/2  and 
1/4  and  1/8  (and  so  on)  the  size  of  the  basic  wavelet.  Whenever 
the  "wavelet  length"  shrinks  by  a  factor  of  two,  the  "wavelet 
frequency"  can  be  thought  of  as  doubling.  This  shrinkage  factor 
is  sometimes  called  scale  or  scale  level.  Wavelet  basis  functions 
are  all  related  by  multiples  of  the  constant  ratio  (2  :  1) .  The 
basic  wavelet  is  computed  from  the  scaling  function.  The 
selection  of  the  scaling  function  determines  all  of  the  remaining 
basis  functions.  A  remarkable  fact  is  that  there  are  an  infinite 
number  of  scaling  functions,  each  of  which  defines  a  complete 
wavelet  basis.  This  provides  tremendous  flexibility  in  selecting 
basis  functions  which  are  appropriate  for  different  systems. 

In  general,  the  computational  complexity  of  wavelet  methods 
is  0(n),  which  means  that  the  computational  complexity  is  of  order 
n.  The  efficiency  is  the  direct  result  of  the  simplicity  of  the 
wavelet  transform  process.  The  process  starts  by  separating  the 
signal  information  in  the  smallest  wavelets  from  information  in 
all  the  larger  wavelets  scales.  Details  are  given  in  Section  2.4. 
The  output  of  the  first  stage  is  processed  again  by  the  same 
method  and  is  repeated  for  each  successive  scale.  This  recursive 
structure  reduces  the  amount  of  data  to  be  processed  at  each 
successive  level  by  a  factor  of  two  which  reduces  the 
computational  cost  for  each  successive  transform  level. 

Information  which  varies  rapidly  over  just  a  few  data  points  is 
separated  from  information  which  varies  over  many  points.  The 
procedure  is  stopped  at  the  largest  wavelet  of  interest . 

The  outputs  of  the  wavelet  transform  are  coefficients  which 
represent  the  similarity  of  the.  signal  (as  a  function  of  time)  to 
the  wavelets  of  different  shrinkage  factors  and  times.  The  output 
of  a  wavelet  transform  can  be  plotted  on  a  grid  that  has  time  on 
the  x-axis  (or  t-axis)  and  shrinkage  factor  on  the  y-axis.  This 
grid  is  called  phase  space  and  is  used  to  graphically  display  the 
relationships  between  signal  information  at  different  shrinkage 
factors  and  times. 

It  is  helpful  to  compare  and  contrast  the  wavelet  transform 
with  the  well  known  Fourier  transform.  The  Fourier  transform  is 
also  invertible,  energy-preserving  and  linear.  The  Fourier  basis 
functions  are  orthogonal  and  complete  for  the  entire  set  of  finite 
energy  discrete  signals. 

The  Fourier  transform  separates  signal  information  by 
frequency.  The  Fourier  transform  requires  an  a  priori  choice  of 
input  data  block  size.  The  Fourier  basis  functions  are  constant 
frequency  complex  exponential  functions  each  of  which  persists  as 
long  as  the  block  size.  The  basis  functions  are  uniformly  spaced 
in  frequency.  The  frequencies  are  separated  by  a  constant 
interval  rather  than  a  constant  ratio.  Since  all  of  the  basis 
functions  are  as  long  as  the  input  data  block,  they  all  have  the 
same  (lack  of)  time  resolution.  The  output  of  a  Fourier  transform 
contains  information  about  how  the  energy  in  the  signal  is 
distributed  among  the  frequencies  in  the  signal.  However 
information  about  how  the  energy  is  distributed  in  time,  about 
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when  it  occurred,  is  not  available  in  the  Fourier  transform 
representation.  All  that  can  be  inferred  is  that  the  frequency 
was  present  somewhere  in  the  block  and  what  fraction  of  the  signal 
energy  it  accounted.  The  computational  complexity  of  the  fast 
Fourier  transform  (for  the  commonly  used  Cooley-Tukey  algorithm) 
is  0 (nlog2 (n) ) . 

The  Fourier  transform  separates  signal  information  into 
uniformly  spaced  frequency  components.  The  Fourier  transform  has 
fine  frequency  resolution  and  a  complete  lack  of  time  resolution. 
Increasing  the  input  data  block  size  increases  the  range  of 
frequencies  which  the  Fourier  method  can  resolve,  but  decreases 
the  time  information  available  from  a  signal. 

In  summary,  the  primary  difference  between  the  Fourier 
transform  and  the  wavelet  transform  is  in  how  each  separates 
signal  information  between  time  and  frequency  or  a  frequency 
related  parameter,  shrinkage  factor.  Wavelet  transforms  separate 
signal  information  by  shrinkage  factor  and  time.  The  number  of 
shrinkage  factors  used  and  the  number  of  times  resolved  are 
jointly  limited  by  the  input  data  block  size.  There  are  an 
infinite  number  of  wavelet  basis  from  which  an  appropriate  basis 
can  be  selected.  The  computational  complexity  of  the  wavelet 
transform  is  0(n),  less  than  the  Fourier  transform  0(nlog2(n)). 
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2 . 1  The  Scaling  Function,  Wavelets,  and  the  Wavelet 

Transform 

Wavelet  methods  separate  the  components  of  a  signal  by  time 
and  a  shrinkage  factor  that  is  related  to  frequency.  With  a 
wavelet  transform,  both  the  time  resolution,  or  "correlation 
length",  and  the  shrinkage  factor  resolution  vary  logarithmically. 

The  wavelet  technique  takes  into  account  the  reciprocal 
relationship  between  time  and  frequency  (or  any  type  of  structure 
which  is  expressed  across  multiple  data  points) .  To  identify  or 
locate  the  position  of  a  particular  shape,  such  as  an  oscillation, 
in  a  set  of  data,  one  must  look  for  relationships  among  the  data 
values;  a  structure  or  an  oscillation  exists  only  across  a  set  of 
data,  and  not  in  a  single  point  value.  A  number  has  no  frequency, 
no  structure.  Conversely,  properties  which  exist  throughout  a  set 
of  data  cannot  be  said  to  have  a  particular  location  within  that 
set.  Wavelet  signal  representation  techniques  take  this  trade-off 
into  account  by  allowing  small  sets  of  data  to  be  combined  and 
correlated  to  derive  structural  or  shape  information  about  that 
subset,  without  requiring  a  complete  transformation  of  the  signal 
into  a  particular  type  of  structural  information,  the  way  a 
Fourier  transform  does.  Thus,  unlike  the  Fourier  transform  one  is 
allowed  to  exchange  a  small  amount  of  temporal  resolution  for  a 
small  amount  of  information. 

It  is  no  accident  that  these  properties  are  reflected  in  the 
characteristics  of  many  natural  signals.  Signals  with  time 
varying  characteristics,  like  speech,  music,  seismic  signals  and 
underwater  acoustic  signals  are  all  best  analyzed  by  a  system 
capable  of  resolving  both  frequency  and  time.  Furthermore,  many 
signal-producing  phenomena  have  octave  band  structure  due  to  the 
presence  of  harmonics  within  the  signal  and  respond  well  to 
wavelet  analysis.  Transient  events  also  respond  well  to  wavelet 
analysis  in  that  the  identification  of  precisely  located 
phenomena,  such  as  the  sharp  onset  of  a  signal,  requires  the 
ability  to  resolve  its  location  in  time  with  a  very  short 
wavelength,  while  the  characteristics  of  later,  more  persistent, 
parts  of  the  signal  may  require  the  ability  to  identify  longer 
wave  shapes . 

Scaling  functions,  wavelets,  and  the  wavelet  transform  are 
discussed  in  more  detail  in  the  next  three  sections  and  also  in 
Appendix  B. 
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2 . 2  The  Scaling  Function 

The  scaling  function  is  at  the  core  of  any  wavelet  based 
representation  of  a  signal.  We  will  discuss  only  compactly 
supported  wavelets  in  this  report.  The  scaling  function  has  three 
essential  properties.  The  first  is  that  it  is  compactly 
supported.  This  means  that  the  scaling  function  is  exactly  zero 
outside  a  bounded  region  of  the  real  line.  The  scaling  function 
is  only  locally  non  zero. 

The  second  essential  property  is  that  the  scaling  function  is 
orthogonal  to  integer  translates  of  itself.  The  importance  of 
this  will  become  clear  a  little  later.  The  third  property  is  that 
the  scaling  function  is  intimately  related  to  smaller,  or  scaled 
versions  of  itself.  This  relationship  is  expressed  concisely  by 
the  scaling  equation: 

N-1 

<p(x)  =  ^  ait<p(2x  -  k)  (2.1) 

where  <p(x)  is  the  scaling  function.  The  function  <p(2x)  is  a 
smaller,  scaled  down  (by  a  factor  of  two),  version  of  <p(x)  .  The 
scaling  equation  states  that  <p(x)  is  equal  to  a  weighted  sum  of 
these  small  versions  of  itself.  The  numbers  oi^,  of  which  only 
finitely  many  (N,  which  is  an  even  number  here)  are  non  zero,  are 
called  the  scaling  coefficients.  N  is  the  size  of  the  wavelet 
system.  The  support  of  <p,  the  region  on  which  it  is  non  zero,  is 
the  interval  [0,  N-1]  .  The  coefficients  a)c  must  satisfy  certain 
conditions  in  order  for  the  scaling  function  to  exist  and  satisfy 
the  scaling  equation.  There  turns  out  to  be  an  infinite  number  of 

sets  of  scaling  coefficients  for  every  even  N  >  2.  It  is  the 

choice  of  the  oi^,  from  among  this  set,  which  determines  the 
detailed  shape  of  <p(x)  . 

Scaling  functions  are  commonly  selected  from  the  class  of 
Daubechies  functions,  which  have  several  important 
characteristics.  They  are  relatively  smooth  and  have  certain 
approximation  properties  (i.e.,  vanishing  moments).  These  systems 
will  be  referred  to  as  D2  (which  is  also  the  "Haar”)  D4,  D6, 
D8...Dn,  where  n  is  the  size  of  the  system.  The  first  four 
Daubechies  scaling  functions  are  shown  in  Figure  2.1. 

The  scaling  function  is  the  basic  unit  from  which  a  level  of 

detail  is  constructed.  This  is  done  by  considering  the  set  of 
functions  which  can  be  represented  as  a  linear  combination  of 
shifted  versions  of  the  scaling  function.  That  is,  we  define  a 
collection  of  functions  at  "scale  level"  j,  which  we  write  Vj  to  be 
the  set  of  functions  which  are  linear  combinations  of  the 
functions  «p(2jx  -  k)  ,  where  k  is  an  integer.  The  factor  2i 
multiplying  x  has  the  effect  of  shrinking  the  support  of  tp  to  the 
interval  from  0  to  (N  -  l)/2),  and  the  shift  by  k  moves  these  small 
functions  around.  Thus  a  function's  components  at  scale  level  j 
are  expressed  by  the  equation: 
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f  j  (X) 


(2.2) 


=  X  Cj.ic<l>(2)x  -  k) 
k£Z 

where  fj  is  the  part  of  f  resolvable  at  the  scale  level  j. 

This  idea  can  also  be  expressed  by  stating  that  Vj  is  the 
space  spanned  by  the  set 

{<p(2jx  -  k)  IkeZ)  .  (2.3) 

This  set  of  functions  forms  an  orthonormal  basis  for  Vj .  Functions 
in  Vj  are  uniquely  expressible  as  linear  combinations  of  the  basis 
functions,  and  the  basis  functions  all  have  unit  "energy".  Thus 
the  set  of  functions  {<p(2)x  -  k)  }  form  an  orthogonal  set  of 
"templates"  for  Vj. 

The  effect  of  performing  a  transform  with  such  a  set  of  basis 
functions  is  to  identify,  within  the  signal,  those  parts  or 
components  which  are  similar  to  the  basis  functions  at  the  given 
scale  level.  Similarly,  a  Fourier  transform  has  oscillatory 
functions  as  a  basis,  and  identifies  the  relative  contribution  of 
each  frequency  to  the  overall  signal. 

With  shifted  versions  of  the  scaling  function  as  a  basis,  the 
wavelet  transform  will  identify  components  which  are  similar  to  a 
particular  shifted  copy  of  the  scaling  function;  that  is, 
representations  of  a  function  in  the  scaling  function  basis 
identifies  features  locally  in  time,  since  the  scaling  function  is 
compactly  supported,  and  locally  in  scale  level,  because  the 
scaling  function  has  structure. 

There  is  another,  equally  good,  set  of  orthogonal  "templates" 
for  this  scale  level  Vj.  This  set  of  "templates"  gives  a  different 
type  of  information  than  the  one  presented  above.  In  our  previous 
basis,  all  of  the  resolution  within  the  scale  level  j  was  in  the 
temporal  domain;  each  coordinate  corresponded  to  a  position  in 
time.  The  new  basis  will  trade  some  of  this  temporal  resolution 
for  some  additional  structural  information;  for  each  scale  level 
it  will  give  us  two  sets  of  coefficients,  one  set  which  represents 
large  structure,  while  the  other  set  represents  small  ("fine") 
structure.  The  process  extracts  or  filters  out  the  components  of 
the  scale  level  j  which  cannot  be  regarded  as  part  of  the  coarser 
scale  level,  j  -  1, 
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2.3  Wavelets 


This  idea  of  dividing  the  scale  level  Vj  into  a  coarser  version  of 
itself,  Vj_i,  and  a  difference  space,  which  we  will  call  Wj-i,  can 
be  compactly  expressed  as  an  orthogonal  splitting  of  the  space  Vj 
into  two  perpendicular  spaces  Vj-i  and  Wj-i: 

Vj  =  Vj_i  e  Wj.i.  (2.4) 

The  reason  this  can  be  done  efficiently  comes  from  the  scaling 
equation.  Since  <p(x)  can  be  expressed  as  a  linear  combination  of 
translated  versions  of  «p(2x)  the  coarser  scale  level  Vj_i  is 
contained  within  the  finer  scale  level  Vj: 

Vj-i  c  Vj  (2.5) 

Repetition  of  the  argument  shows  that 

Vj.i  C  Vj  C  Vj^i  C  Vj^2  ...  (2.6) 

The  difference  space,  Wj_i,  contains  all  that  remains  when  the 
coarser  scale  information  is  removed.  However,  since  Wj_i  is 
contained  in  Vj,  it  is  also  expressible  as  a  linear  combination  of 
translates  of  <p(2jx)  . 

The  actual  set  of  functions  that  are  used  to  span  the  space 
Wj-i  are  the  orthonormal  basis  formed  by  the  functions 


V(2jx)  = 

N-1 

^  (-l)KaN-)c-i<p(23+ix  - 

k) 

(2.7) 

or  in  the  case  of  j  =  0, 

V  (x) 

N-1 

=  y  (-l)!«aN-K-i<p(2x  - 

k=0 

k) 

(2.8) 

Notice  that  the  signs  now  alternate  in  the  sum,  and  the  order  of 
the  coefficients  oj^  has  been  reversed  (k  -►  N  -  )c  -  1)  .  These 
changes  make  v(x)  orthogonal  to  «p  (x)  .  The  full  basis  for  Wj  is 
formed  by  taking  shifted  versions  of v,  i.e..  Basis  (Wj)  = 

{2j/2y (2jx-k)  |k  an  integer}.  The  support  of  v(x)  is  easily  seen  to 
be  the  same  as  the  support  of  ^(x),  and  this  is  true  for  the 
shrunken  versions  as  well,  that  is,  the  support  of  v(2jx-k)  is  the 
same  as  the  support  of  v(2jx-k)  .  The  normalization  term  2j/2 
maintains  unit  energy  in  the  functions.  Figure  2.2  shows  the 
wavelets  which  correspond  to  the  scaling  functions  presented  in 
Figure  2.1. 

Thus  the  transformation  from  the  first  representation  of  Vj, 
where  f j  was  expressed  as  a  linear  combination  of  shifted  versions 
of  (p(2)x),  to  the  new  representation,  where  fj  is  expressed  as  a 
sum  of  translates  of  «p(2j-ix)  and  translates  of  v(23-ix),  gives  us 
new  information  about  the  shapes  and  structures,  perhaps 
frequencies,  present  in  fj.  This  is  at  the  expense  of  some  temporal 
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resolution,  because  the  new  basis  functions  are  twice  as  long.  The 
new  basis  incorporates  the  inter-relationships  among  larger 
subsets  of  the  data,  providing  correlative  information.  As  a 
result  of  the  spectral  refinement,  there  has  been  a  loss  of 
temporal  resolution. 

The  scaling  function,  being  the  origin  for  all  of  the  spaces 
Vj,  forms  the  connection  between  these  spaces  via  the  scaling 
equation,  equation  2.1.  The  basic  wavelet,  v,  represents  the 
differences  between  the  scale  levels. 
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2 . 4  The  Wavelet  Transform 


The  exchange  of  temporal  resolution  in  Vj  for  scale  resolution  in 
the  division  of  Vj  into  Vj_i  and  Wj_i  forms  the  basic  unit  of  the 
Wavelet  transform.  Since  the  definition  is  independent  of  scale 
level,  it  can  be  repeatedly  applied  in  the  same  way.  Furthermore, 
the  operations  involved  in  the  transformation  require  only  the 
expansion  coefficients  of  the  function  f (x)  in  the  basis  at  the 
current  scale,  and  not  the  actual  values  of  the  function.  The 
computation  is  very  simple  and  efficient  because  of  the  close  link 
between  the  functions  <p  and  v. 

The  basic  operation  involved  in  a  wavelet  transform  is  the 
conversion  of  temporal  resolution  into  structural  or  spectral 
information.  This  basic  single  step,  essentially  a  filter, 
exchanges  half  the  temporal  resolution  of  a  signal  for  twice  the 
"frequency”  resolution;  the  product  of  the  two  remains  the  same. 
More  importantly,  this  operation  can  be  repeated  to  gain  any 
desired  level  of  detail  in  the  structural  or  spectral  realm,  while 
only  imposing  a  reciprocal  loss  of  resolution  in  the  temporal 
domain.  This  is  in  sharp  contrast  to  Fourier  transform 
techniques,  where  one  either  gets  all  the  available  frequency 
information,  or  none  of  it,  with  no  intermediate  stages  of 
knowledge  available. 

The  wavelet  transform  allows  one  to  move  gradually  between 
the  two  extremes  present  in  the  Fourier  transform,  successively 
gaining  shape  or  structure  information.  In  a  sense,  the  wavelet 
transform  interpolates  between  the  frequency,  or  structure,  domain 
and  the  temporal  domain.  This  step  by  step  transformation  can  be 
understood  in  terms  of  trade-offs  between  relative  time  resolution 
and  relative  frequency  resolution. 

Since  the  transform  is  defined  in  terms  of  operations  on  the 
coefficients  of  the  representation,  and  not  the  actual  values  of 
the  scaling  or  wavelet  functions,  the  output  from  a  single  stage 
of  the  transform  is  exactly  what  the  next  stage  requires  for 
input.  This  easily  pipelined,  recursive  structure  is  what  makes 
the  wavelet  transform  rapidly  computable.  While  many  such 
structures  are  made  possible  by  the  wavelet  transform,  one  in 
particular,  the  one-sided  or  Mallat  transform,  has  proven  to  be 
exceptionally  useful  in  analyzing  signals.  Figure  2.3  shows  one 
operation  of  the  wavelet  transform.  The  second  operation  uses  the 
"low  pass"  (see  Figure  2.4)  output  coefficients  from  the  first 
stage  as  the  input  for  the  second  stage,  etc.  The  number  of  data 
points  at  each  level  is  reduced  by  a  factor  of  two  by  the 
conversion  of  temporal  information  into  spectral  information. 
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Figure  2.3: 


one  Operation  of  the  Wavelet  Transform 


Figure  2.4:  Tree  Structure  for  the  Partial  Sub-band  TVee  Analysis. 
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2 . 5  The  Computational  Complexity  of  Wavelet  Transforms 

This  section  is  concerned  with  the  computational  efficiency  of 
wavelet-based  analysis  techniques.  The  computational  complexity 
of  several  types  of  wavelet  transforms  is  developed  and 
comparisons  are  made  to  the  computational  complexity  of  the  FFT. 
Both  pre-  and  post-processing  requirements  are  ignored.  It  is 
assumed  that  the  output  from  the  wavelet  transform  or  FFT  is  the 
desired  result.  The  input  block  size  is  the  factor  which 
determines  the  computational  cost.  The  first  case  is  a  wavelet 
transform  which  resolves  only  a  portion  R=1/2J  of  the  total 
bandwidth,  and  resolves  it  as  finely  as  possible  using  recursive 
wavelet  techniques.  The  simplest  example  of  such  a  transform  is 
the  familar  Mallat  Transform,  which  calculates  the  decomposition 
of  a  signal  on  the  basis  of  scale,  i.e.,  it  "homes  in"  on  low 
frequencies.  This  fundamental  structure  requires 

OPS(Mallat)  =  a(K)N(l-l/2J)  #  of  operations  (2.9) 

where  a(K)=(K+l)  multiplies  and  (K)  additions  per  output  point.  The 
number  of  input  data  is  N,  and  J  is  the  finest  level  calculated; 
that  is,  the  basic  wavelet  decomposition  operation  is  applied  J 
times.  The  number  of  nonzero  coefficients  is  K,  which  we  have 
also  called  the  length  of  the  wavelet  coefficient  matrix.  Note 
that  whatever  the  depth  of  the  decomposition,  the  operation  count 
never  exceeds  a(K)N.  This  operation  count  also  applies  to  any 
wavelet  transform  which  "zooms  in"  on  a  single  location  in  phase 
space,  allowing  other  side-bands  to  remain  unchanged.  These  are 
not  partial  wavelet  transforms:  each  is  a  complete  representation 
in  a  wavelet  basis.  Each  is,  however,  a  partial  frequency 
decomposition.  That  is  a  powerful  advantage  because  only  what  is 
required  need  be  calculated. 

In  applications  where  one  wishes  to  resolve  frequencies  (or 
some  other  structure)  to  some  pre-specif ied  resolution  (say  1/2'^)  , 
and  one  is  interested  in  a  small  number  of  sub-bands,  wavelets  are 
very  computationally,  efficient .  This  case  includes  the  Mallat 
Transform,  and  any  other  wavelet  processing  scheme  which  generates 
only  a  small  subset  of  the  finest  resolution  cells. 

If  the  application  requires  that  the  temporal  resolution 
existing  in  a  subband  be  converted  to  frequency  or  structure 
information,  the  wavelet  transform  can  be  further  applied  in  a 
uniform  fashion  to  derive  this  information,  as  in  Figure  2.4. 

Since  the  formula  for  this  case,  equation  2.10,  is  rather 
complicated,  comparisons  were  made  with  the  complexity  of  an  FFT 
calculation  for  the  full  bandwidth  of  the  signal,  and  for 
consistency  the  wavelet  transform  was  also  carried  to  its  full 
resolution  within  the  resolved  subband.  Since  there  are 
essentially  no  savings  possible  with  an  FFT  for  resolving  only  a 
portion  of  all  subbands,  the  full  complexity  of  3n+nlog2n  was  used 
for  comparison. 
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If  we  define 


1.  N  =  the  number  of  data  points; 

2 .  K  =  the  size  of  the  wavelet  filter  (assumed  to  be 
constant  throughout  the  procedure)  and  a(K)  is  the 
computational  cost  for  the  fixed  filter  length; 

3.  R  =  the  portion  of  the  bandwidth  of  the  signal  which  is 
resolved  (assumed  to  be  1/2^  for  some  L) ; 

4.  J  =  the  finest  level  of  resolution  (the  width  of  the 
frequency  bands  resolved  is  1/2^) , 

then  the  number  of  operations  required  to  generate  the  desired 
results,  along  with  all  of  the  remai.iing  results  that  are  required 
to  have  a  complete  representation,  is: 

a(K)N(l+R(J-2+log2 (R)  ) /2)  (2.10) 

The  comparison  is  summarized  by  Figures  2.5  and  2.6,  (o(K)  counts 

the  number  of  real  multiplies)  which  show  the  filter  length  for 
which  the  two  techniques  are  ecpaally  complex  for  a  variety  of 
values  of  R.  See  Figure  2.4  for  an  example  of  such  a 
decomposition . 

Figures  2.5  and  2.6  relate  increased  wavelet  computationals 
due  to  increased  number  of  coefficients,  K,  to  wavelet  superiority 
versus  FFT  (in  terms  of  computations  required),  which  increases  as 
the  data  bloclc  size  increases. 
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Figure  2.5:  Computational  complexity  of  the  Wavelet  Transform  vs.  the 
FFT:  Filter  Length  vs.  N  for  equivalent  complexity  for  various  values  of  R. 


32  bins.  R  =  32/N. 
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3.0  Basic  Radar  System  Operation 

Radar  is  an  active  sensor  system  that  transmits  electro¬ 
magnetic  energy  and  interprets  the  echoes  reflected  from 
objects.  Objects  in  the  propagation  path  scatter  this  energy 
in  directions  determined  by  the  physical  and  electromagnetic 
characteristics  of  the  object.  The  reflected  energy  is  then 
received,  processed,  and  perhaps  interpreted  by  the  radar 
system.  Since  electromagnetic  energy  propagates  with  a 
(nearly)  constant  velocity  in  the  atmosphere,  the  round  trip 
travel  time  t  from  the  radar  to  an  object  and  back  again  can 
be  converted  to  distance  by  the  formula 

r  =  ct/2  (3.1) 

where  r  is  the  range  to  the  object  and  c  is  the  speed  of 
1 i ght . 

Radars  are  built  for  many  purposes.  Some  are  used  to 
detect  objects  such  as  aircraft,  ships,  or  tanks.  Others  are 
used  to  provide  weapon  guidance  information.  Still  others 
are  used  to  map  the  surface  of  the  earth  or  for  weather 
forecasting.  In  the  first  two  cases,  users  are  interested  in 
the  energy  reflected  from  the  "target"  and  they  desire  that 
it  be  separated  from  energy  received  from  the  earth  or 
sources  of  clutter.  Similarly,  a  "target"  radar  signal  may 
be  corrupted  by  "noise"  which  is  undesired  or  extraneous 
energy  received  with  the  sighal  of  interest.  The  total  noise 
level  present  in  a  radar  depends  on  many  factors  including 
the  weather,  operation  of  nearby  electrical  equipment  (either 
friendly  or  hostile)  and  many  other  factors.  The  radar  data 
processing  requirements  will  vary  from  system  to  system 
depending  on  the  intended  use. 

Radar  processors  may  have  to  perform  any  or  all  of  the 
following  tasks.  "Detection"  means  that  the  radar  has  found 
some  object  that  may  possibly  be  a  target  based  upon  some 
very  fast  data  processing  procedure  designed  to  eliminate 
most  noise  and  homogeneous  clutter.  One  such  process  is 
called  CFAR  (for  "constant  false  alarm  rate")  where  thresh¬ 
olds  are  set  to  eliminate  desired  percentages  of  noise. 
"Discrimination"  is  the  process  of  separating  detected 
targets  from  detected  clutter.  "Classification"  is  the 
process  of  identifying  particular  types  of  targets  (some  of 
which  may  be  more  important  as  a  missile  target  than  others  - 
for  example,  tanks  versus  trucks).  Also,  the  radar  may 
perform  "tracking"  which  includes  computing  and  frequently 
updating  target  position  and  motion  information. 
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The  Millimeter  Wave  Radar  Data  Base 


3.  1 


Millimeter  wave  (MMW)  radar  systems  are  particularly 
well  suited  for  use  in  tactical  missile  systems  because  they 
are  physically  small  and  their  short  wavelengths  provide 
excellent  target  range  resolution  in  nearly  all  weather 
conditions.  Missile  sensors  are  usually  used  only  once;  so, 
they  must  be  inexpensive.  Also,  MMW  radars  have  some  resist¬ 
ance  to  electronic  countermeasures  (i.e.,  jamming).  A 
typical  MMW  radar  block  diagram  is  given  in  Figure  3.1, 
and  areas  that  wavelets  methods  might  apply  to  are  high¬ 
lighted  in  Figure  3.2. 

The  radar  data  used  in  studies  on  this  contract  was 
taken  by  a  35  Ghz  MMW  radar.  This  data  base  contained 
coherent  radar  looks  at  a  trihedral  (which  is  visible  when 
the  transmitter  -  receiver  polarization  is  "odd"),  looks  at 
two  dihedrals  (which  are  visible  when  the  transmitter  - 
receiver  polarization  is  "even")  separated  in  range,  looks  at 
many  360  degrees  of  azimuth  aspect  angles  of  an  M60  tank,  and 
an  M35  truck.  The  trihedral  (representing  odd  bounce 
scattering)  is  visible  when  the  transmitter-receiver  are 
"odd"  or  of  opposite  sense  (transmitting  right  hand  circular 
polarized  energy  and  receiving  lef*^  hand  circular  polarized 
energy  or  vice-versa).  The  dihedrals  (representing  even 
bounce  scattering)  are  visible  when  the  transmitter-receiver 
polarizations  are  "even"  or  of  the  same  sense  (transmitting 
right  and  receiving  right  or  transmitting  left  and  receiving 
left) . 

Each  polarization  of  a  radar  look  consists  of  63 
frequency  stepped  pulses  transmitted  at  and  received  from  the 
same  target  or  area  of  ground.  A  64  point  (one  zero  added) 
FFT  then  forms  a  matched  filter  with  the  received  radar 
frequency  stepped  pulse  train.  That  is,  phase  change  (due  to 
a  point  scatterer)  in  the  pulse  train  will  "match"  one  of  the 
columns  of  a  discrete  Fourier  transform  ( DFT )  matrix.  In 
that  way,  the  original  pulse  length  that  was  twice  63  feet 
(for  round  trip  travel)  is  subdivided  to  about  one  foot  range 
resolution  by  FFT  processing.  The  result  is  called  a  high 
range  resolution  profile  (HRRP).  This  method  of  achieving 
high  range  resolution  allows  the  use  of  fairly  long  pulses. 
That  in  turn  allows  long  range  target  detection  with  low  peak 
power  radar  transmitters.  There  are  two  "odd"  polarization 
HRRPs  and  two  "even"  polarization  HRRPs  available  from  each 
radar  look.  Additional  details  of  this  data  base  are  given 
in  Appendix  D. 

Typical  (head-on  aspect)  HRRPs  and  partial  wavelet 
reconstructions  using  the  D2  (Haar)  wavelet  basis  and  50%, 
25%,  and  10%  of  the  coefficients  to  reconstruct  the  original 
(100%)  HRRP  are  shown  in  Figures  3.3  through  3.6.  Figure  3.3 
is  a  plot  of  the  results  for  an  even  polarization  M60  tank. 
Figure  3.4  is  a  plot  for  the  odd  polarization  M60.  Figure 
3.5  is  a  plot  of  the  results  for  an  even  polarization  M35 
truck.  Figure  3.6  is  a  plot  of  the  results  for  an  odd 
polarization  M35.  It  can  be  seen  from  these  figures  that 
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reasonably  fine  structure  is  preserved  along  with  the  target 
extent,  even  when  only  25%  of  the  wavelet  coefficients  are 
used  to  reconstruct  the  original  target  image.  This  reduces 
target  image  dimension  by  a  factor  of  four.  It  is 
conjectured  that  features  from  the  reconstructed  image  may  be 
more  robust  than  features  from  the  original  image.  The  next 
phase  of  this  contracc  will  investigate  this  conjecture  in 
detai 1  . 
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MMW  RADAR  BLOCK  DIAGRAM 


WAVELET  IMPACT  AREAS 


Figure  3.4:  D2  RECONSTRUCTED  M60  -  ODD 
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3 . 2  Wavelet  Approximation  to  the  FFT 

Several  studies  were  conducted  using  the  radar  data  to 
develop  and  demonstrate  wavelet  processing  methods  for 
approximating  the  Fourier  transform  and  extracting  target 
features.  The  experiments  ranged  from  simple  demonstrations  of 
the  methods,  which  proved  that  they  could  be  performed,  to  a 
series  of  trials  to  evaluate  the  sensitivity  of  the  methods  to  the 
choice  of  wavelet  basis  functions  and  target  signature  variations. 
The  sensitivity  of  these  methods  to  the  selection  of  basis 
functions  is  important,  in  part  because  the  choice  of  basis 
function  directly  influences  the  computational  complexity  of  the 
method.  The  sensitivity  to  signature  variation  is  important 
because  the  military  targets  (i.e.,  tank  and  truck)  exhibit 
significant  signature  fluctuations  as  a  function  of  viewing  angle. 

This  section  reports  the  results  of  these  studies  which  were 
performed  to  assess  the  basic  feasibility  of  performing  wavelet 
approximations  to  the  discrete  Fourier  transform  (DFT) .  The 
mathematical  details  of  wavelet  transform  theory  are  discussed  in 
Appendix  B.l,  and  details  of  wave  transform  approximations  to 
Fourier  transforms  are  in  Appendix  B.2. 

The  method  we  chose  to  quantify  the  error  of  approximation  is 
the  ratio  of  the  mean  squared  error,  to  the  total  energy  in  the 
signal.  This  method  is  sensitive  to  any  loss  of  target  energy  by 
the  approximation  method  and  the  normalization  allows  the  direct 
comparison  of  different  signals. 

Figure  3.7  illustrates  the  normalized  error  of  approximation 
as  a  function  of  the  number  of  wavelet  terms  retained  in  the 
wavelet  FFT  approximation  formula  (B.30) .  Note  that,  for  this  M60 
tank  scattering  profile,  the  error  of  approximation  decreases  very 
rapidly  in  the  first  5  or  6  terms  followed  by  a  more  gradual  and 
steady  convergence  to  zero.  This  behavior  is  typical  of  the 
method  applied  to  other  target  profiles,  as  illustrated  in  Figure 
3.8.  The  convergence  of  this  method  for  each  of  the  target 
profiles  (which  vary  in  aspect  or  viewing  angle)  is  rapid  and 
consistent.  The  method  shows  little  sensitivity  to  changes  in 
aspect  angle  and  therefore  to  target  signature  variations. 

Figure  3.9  is  a  plot  of  the  error  for  a  number  of  wavelet 
bases.  The  functions  considered  ranged  from  a  2  coefficient 
wavelet  basis  (the  Haar  basis)  to  an  eight  coefficient  wavelet 
basis  (Daubechies  -  8.)  There  is  very  little  difference  in  their 
performance.  This  surprising  result  indicates  that  the 
approximating  wavelet  basis  does  not  have  to  be  chosen  very 
carefully  and  that  the  computational  advantages  of  short  wavelets 
can  be  fully  exploited  in  this  application. 

In  the  previous  cases,  the  wavelet  terms  which  were  retained 
in  the  approximation  were  the  wavelet  coefficients  with  the 
largest  absolute  value.  They  were  determined  by  sorting  the 
wavelet  coefficients.  Sorting  is  an  expensive  computational 
process  which  would  quickly  eliminate  any  computational  advantage 
wavelets  would  have  relative  to  Fourier  methods.  There  is  an 
alternative  formulation  of  the  wavelet  approximation  to  the 
Fourier  transform  which  uses  all  of  the  wavelet  coefficients  at  a 
given  scale  (Eq.  B.33)  rather  than  coefficients  selected  by 
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Figure  3.7;  Error  of  Approximation  versus  Number  of  Wavelet  Terms. 


Figure  3.8:  Error  of  Approximation  versus  Number  of  Wavelet  (0)  Terms  for 
Selected  Target  Profiles. 


Number  of  \|/ Terms 

Figure  3.9:  i^rror  of  Approximation  versus  Number  of  Wavelet  (rjf)  Terms  for 
Selected  Wavel^'t  Bases. 


Number  of  \|/  Levels  (J's) 

Figure  3.10:  Error  of  Approximation  versus  Scale  (0)  Levels  for  Selected 
Wavelet  Bases. 
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magnitude .  An  experiment  was  performed  to  evaluate  the 
feasibility  of  eliminating  entire  scales  from  the  approximation. 
Figure  3.10  shows  the  results  for  several  wavelet  basis  functions. 
The  approximation  is  very  poor  until  the  smallest  scale  wavelet 
coefficients  are  included.  This  means  that  the  energy  was 
distributed  in  wavelet  coefficients  across  all  scales. 

The  results  of  the  previous  experiment  motivated  a  search  for 
a  method  to  preprocess  the  radar  data  so  that  the  target  energy 
could  be  placed  into  predictable  scale  levels  and  thereby 
eliminate  the  need  to  select  (or  sort)  the  most  significant 
wavelet  coefficients.  This  led  us  to  adapt  a  method  called  "phase 
unwrapping"  to  this  problem.  Appendix  C  describes  the 
mathematical  details  of  the  procedure.  When  it  is  applied  to  the 
radar  data,  the  target  energy  is  concentrated  into  the  lowest 
scale  levels  and  the  approximation  method  can  discard  the  highest 
scale  level  terms  with  little  loss  of  information.  This 
significantly  reduces  the  computational  complexity  of  the  method. 
Figure  3.11  illustrates  the  reduction  in  error  it  produced. 

In  conclusion,  it  is  feasible  to  apply  the  wavelet 
approximation  to  the  Fourier  transform  to  the  millimeter  wave 
sensor  problem  of  producing  HRRPs.  The  error  of  approximation  can 
be  controlled  to  any  level  desired  by  adjusting  the  number  of 
wavelet  terms  retained  in  the  approximation.  The  computational 
complexity  of  the  wavelet  methods  for  wavelet  support  lengths  of  6 
or  less  are  comparable  to  the  Fast  Fourier  Transform  for  the  radar 
data  used  in  this  study,  which  is  processed  in  blocks  of  63 
complex  data  points.  The  target  data  is  contained  in  a  relatively 
large  fraction  (1/2  to  1/3)  of  the  output  data  points  which 
requires  that  a  large  fraction  of  the  wavelet  coefficients  be 
retained  in  the  approximation.  This  increases  the  computation 
required  for  the  wavelet  methods.  The  target  extent  varies  from 
about  10  fine  range  resolution  cells  to  30  range  resolution  cells. 
The  error  of  approximation  was  about  0.10  (Norm.  MSE)  when  16 
terms  were  retained  and  was  less  than  0.05  (Norm.  MSE)  when  32 
terms  were  retained.  It  should  be  noted  that  the  8  term 
approximations  contain  over  80%  of  the  original  target  energy,  but 
have  relatively  severe  "shape"  distortion. 

Since  the  wavelet  methods  offer  no  computational  advantage 
and  produce  some  error,  they  are  probably  not  an  appropriate 
choice  for  this  particular  radar  data  processing  problem  due  to 
the  small  number  of  data  points  to  process.  However,  if  the  data 
came  from,  for  example,  an  FMCW  (with  a  long  pulse  length)  that 
performed  1024  point  FFT's,  yielding  one  foot  range  resolution, 
tank  and  truck  detections,  which  are  10  to  30  feet  long,  depending 
on  aspect,  could  be  very  efficiently  processed  by  wavelet  FFT 
approximation.  Additional  comparison  results  are  given  in 
Appendix  B.3. 
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Normalized  Mean  Squared  Error 


Number  of  \j/ Levels  (J's) 


M  gure  3.11;  Error  of  Appronmation  versus  Scale  (0)  Leveb  for  Phase  Un¬ 
wrapped  Target 
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3 . 3  Target  Length  Estimation 

This  section  presents  the  results  of  experiments  conducted  on 
the  radar  data  to  develop  and  investigate  the  capabilities  of 
wavelet  based  feature  extraction  methods.  Of  particular  interest 
was  the  feasibility  of  performing  feature  extraction  with  a 
wavelet  based  approach  and  the  stability  of  the  method  with 
respect  to  target  signature  variations  due  to  change  of  aspect 
angle . 

The  feature  which  was  selected  for  analysis  was  target  extent 
which  is  the  length  of  the  target  image  in  the  high  resolution 
radar  range  data.  This  feature  is  important  because  it  requires 
that  both  the  "front"  and  "back"  of  the  target  be  determined. 
Target  classification  algorithms  need  this  information  to  delimit 
the  region  in  the  high  resolution  radar  data  which  contains  the 
target  detail.  The  performance  of  the  wavelet  methods  is  compared 
to  a  baseline  method  described  in  appendix  E. 

The  following  procedure  was  used  to  measure  the  stability  of 
the  target  extent  methods. 

1 .  Calculate  the  target  extent  for  a  large  number  of  aspect 
angles.  The  entire  target  data  sets  available  to  the 
subcontractor  were  used  for  the  results  presented  in  this 
report . 

2.  Divide  the  results  into  contiguous  blocks  of  50  items  and 
compute  the  average  extent  for  each  of  the  blocks.  This 
is  interpreted  as  the  true  target  extent  in  this 
neighborhood. 

3.  Calculate  the  differences  between  the  extents  calculated 
in  step  1  and  the  local  averages  calculated  in  step  2 . 
From  the  differences,  calculate  the  standard  deviation 
for  each  block  of  50  estimates.  This  estimate  of  the 
standard  deviation  (locally)  is  a  measure  of  the 
stability  of  the  method. 

4 .  Calculate  the  mean  and  standard  deviation  of  the  local 
standard  deviations  calculated  in  step  3  for  the  entire 
available  target  data  set. 

The  wavelet  algorithm  used  for  extracting  the  target  extent 
was  based  on  thresholding  the  squared  modulus  of  the  large  scale 
component  of  the  one-sided  wavelet  transform  of  the  complex  radar 
image.  This  is  a  method  for  removing  noise  from  the  target  data. 
It  has  two  variations  depending  on  the  order  in  which  mathematical 
operations  are  performed: 

•  IWT(z„)  |2 

•  WT|  (Z;,)  |2 

where  WT ( )  indicates  a  wavelet  transform  and  I  1 2  the  squared 
modulus . 
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The  first  case  requires  the  wavelet  transform  of  complex 
Fourier  transform  data  while  the  second  case  only  requires  that  a 
real  wavelet  transform  be  performed.  The  performance  of  the 
second  method  was  found  to  be  better,  as  well  as  less  costly.  It 
is  described  below. 

The  thresholding  procedure  is  similar  to  the  method  used  in 
the  baseline  method  (described  in  appendix  E)  except  that  only  one 
half  of  the  number  of  high  resolution  range  cells  are  used  to 
compute  the  noise  floor  because  the  one-sided  transform  components 
are  sampled  at  half  the  original  rate,  i.e.  they  are  half  the 
length  of  the  original  sequence.  Two  additional  variations  of  the 
wavelet  transform  method  were  studied  to  investigate  the 
sensitivity  to  scale.  The  one-sided  wavelet  transform  method  was 
recursively  applied  to  the  large  scale  components  resulting  in 
data  which  were  sampled  at  1/4  and  1/8  the  original  rate.  These 
cases  used  only  four  and  two  high  resolution  range  cells, 
respectively,  to  estimate  the  noise  floor.  The  general  wavelet 
method  used  to  estimate  target  extent  is  the  following. 

1.  Compute  the  discrete  Fourier  transform  (DFT)  of  the 
phase  history  data  sequence  z^. 

2.  Compute  the  squared  modulus  of  the  transformed  data  |  Z)<  1 2 . 

3.  Compute  the  large  scale  component  of  the  one-sided 
wavelet  transform  of  |Z(c|2. 

4.  Establish  a  threshold  T>0.  (T  became  25.) 

5.  Compute  the  average  power  A  of  the  8  (or  4  or  2)  high 
resolution  range  cells  furthest  away  from  the  target. 
Since  the  target  is  known  to  be  located  in  the  center  of 
the  range  gate,  these  are  the  first  and  the  last  4  (or  2 
or  1)  high  resolution  range  cells. 

6.  Compute  the  indices  Ii  and  I2  of  the  first  and  last 
entries  of  I  Z)j  1 2  that  exceed  T  x  A 

7.  Compute  target  extent  as  2  x  (I2  -Ii  +  1)  [or  4  x  (I2  -  Ii 
+1)  or  8  X  (I2  -Ii  +1) ] . 

The  computational  cost  in  real  multiplies  of  these  three 
filtering  methods  are  (where  K  is  the  size  of  the  wavelet  system) : 

1.  For  the  one  half  scale  method : Cost (N)  =  (K+1) (N/2)  =* 

0(N)  ; 

2.  For  the  one  quarter  scale  method:Cost (N) = (K+1) (3N/4)  => 
0(N)  ; 

3.  For  the  one  eighth  scale  method:Cost (N) = (K+1) (7N/8) 

=»0(N)  . 
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Section  2.5  contains  a  detailed  discussion  of  computational 
complexity . 

Figures  3.12  through  3.15  are  plots  of  mean  target  extent 
deviation  versus  scale  for  a  threshold  value  of  25.  The  value 
plotted  at  the  64  cell  resolution  position  represents  the  baseline 
method.  The  points  plotted  at  the  32,  16  and  8  resolution 
positions  represent  the  wavelet  methods  at  1/2,  1/4  and  1/8  scale 
respectively.  Note  that  the  wavelet  methods  are  about  as  stable 
as  the  baseline  method  with  resolution  (dimensionality)  reduced  by 
a  factor  of  4 . 
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Resolution  (•  of  Cells) 


=  5  gure  3 .  12:  Mean  Deviation  of  M35  Truck  Target  Extent  -  Wavelets  8,16,  32 
and  Baseline  64  cells  -  LL  Polarization. 


Mean  Deviation  of  Target  Extent,  M35  Truck,  Pol  •  LR 


Resolution  (♦  of  Cells) 

figure  3,13:  Mean  Deviation  of  M35  Truck  Target  Extent  -  Wavelets  8,16, 
32  and  Baseline  64  cells  -  LR  Polarization. 


Mean  Deviation  of  Target  Extent,  M60  Tank,  Pol  -  LL 


-igure  3.14:  Mean  Deviation  of  M60  Tank  Target  Extent  -  Wavelets  8,16,  32 
and  Baseline  64  cells  •  LL  Polarization. 


Mean  Deviation  of  Target  Extent,  M60  Tank,  Pol  «  LR 


Resolution  (•  of  Cells) 

- 1  gu re  3 . 1 5 :  Mean  Deviation  of  M60  Tank  Target  Extent  -  Wavelets  8,16,  32 
and  Baseline  64  cells  -  LR  Polarization. 
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3.4  Wavelet  Target  Classification  Results 

A  two  class  target  classification  experiment  was 
conducted  to  see  if  it  is  feasible  to  improve  radar  target 
classification  signal  processing  by  wavelet  methods,  either 
by  improving  classifier  percent  correct  classification  ( PCC ) 
in  the  presence  of  noise  (added  to  the  very  clean  radar 
target  data)  or  by  decreasing  the  number  of  computations 
required  by  the  classifier  processor,  or  both.  The 
classification  problem  examined  was  for  two  vehicles,  an  M60 
tank  and  an  M35  truck.  The  tank  might  be  a  target  for  a 
missile,  but  the  truck  might  be  a  target-like  object  that  one 
would  not  want  to  waste  a  missile  on.  Classifier  overall  PCC 
was  improved  in  realistic  si gnal / i nterference  situations  by 
up  to  5  percentage  points,  and  computations  were  decreased  by 
a  factor  of  10. 

Quadratic  classifiers  were  used  with  either  odd  degree 
(1,3,5,...)  aspects  as  the  "training”  set  and  even  degree 
(0,2,4,...)  aspects  as  the  "test"  set  or  the  reverse.  Noise 
was  never  added  to  the  "training"  set.  Because  this  is 
actual  radar  data  taken  on  an  outdoor  range  with  very  light 
clutter  background,  the  "clean"  data  has  a  signal  to  clutter 
ratio  of  greater  than  20dB.  The  wavelet  basis  selected  was 
the  Haar  or  Daubechies  2  wavelet  basis.  The  features  used  by 
the  classifier  were  particular  wavelet  coefficients  produced 
by  the  target  HRRPs.  Wavelet  coefficients  were  selected 
logically  but  by  no  means  optimally. 

Specific  results  are  given  in  Figure  3.16.  These 
results  clearly  indicate  that  it  is  feasible  to  improve  radar 
classifier  processors  by  wavelet  methods. 

This  study  used  one  odd  polarization  and  one  even 
polarization.  This  would  be  the  case  if  one  transmitted  only 
one  of  the  two  circular  polarizations,  but  received  both 
circular  polarizations  or  the  reverse.  The  starting  point 
for  this  study  was  the  target  HRRPs,  which  contained  one  foot 
range  resolution  target  magnitude  data.  It  is  assumed  that 
target  length  (as  seen  by  the  radar)  is  no  more  than  32  feet 
In  reality,  target  length  varies  considerably  depending  on 
target  aspect  angle.  For  example,  targets  are  longest  when 
viewed  from  a  head-on  (0  degree  aspect)  or  from  a  180  degree 
aspect . 

In  Figure  3.16 
S  =  total  signal  power 
N  =  total  power  of  noise  added 

S/N  given  is  for  the  tank  (higher  for  the  truck) 

HR54  =  32  odd  &  32  even  HRRP  components 

W64-11  =  11  coefficients  from  HR64  wavelet  transform. 

As  part  of  the  next  phase,  it  will  be  important  to 
perform  a  three  class  classifier  experiment.  This  was  not 
accomplished  during  this  phase  due  to  cost  and  schedule 
constrai nts . 
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Figure  3.16;  Target  Classification  -  Original  vs.  Wavelet 
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3.5  Fine  Range  Resolution  by  Wavelets 

It  is  theoretically  possible  to  obtain  fine  range 
resolution  HRRPs  by  using  wavelet  transforms  instead  of  FFTs 
if  the  waveform  transmitted  is  modified  to  be  appropriate  for 
wavelet  transforms.  This  possibility  is  important  because 
a  radar’s  digital  signal  processor  (DSP)  is  typically  one  of 
the  more  costly  parts  of  the  radar  due  to  heavy  processing 
load  requirements  -  primarily  due  to  the  need  to  perform  many 
FFTs.  Since  a  wavelet  transform  is  considerably  faster  than 
an  FFT,  significant  DSP  cost,  weight,  and  size  savings  will 
result  from  wavelet,  rather  than  FFT,  HRRPs. 

A  wavelet  transform  matrix  W  can  be  constructed  from  the 
product  of  a  number  of  wavelet  "butterfly"  matrices.  This  is 
similar  to  computing  a  discrete  Fourier  transform  (DFT) 
matrix  as  the  product  of  FFT  "butterfly"  matrices.  Unlike 
the  DFT  matrix,  the  wavelet  matrix  W  consists  of  real  (rather 
than  complex)  numbers,  for  the  four  wavelets  of  initial 
interest,  D2  (Haar),  D4 ,  D6 ,  and  D8.  The  column  vectors,  wi , 
W2 ,  ...  ,  wn  of  the  n-by-n  matrix  W  are  orthonormal.  If  y 
denotes  the  column  vector  of  wavelet  coefficients  for  an 
input  column  vector  x,  then  y  could  be  computed  by  (more 
efficiently  by  "butterflies") 

y’  =  x’  W  (3.2) 

and  the  absolute  value  of  the  components  of  y  will  be  a 
fine  range  resolution  HRRP  if  the  n  transmitted  waveforms 
correspond  to  the  wavelets  in  the  following  ways. 

1.  The  n  sequential  transmitted  waveforms  look  like  the  n 
rows  of  the  matrix  W.  (The  orthonormal  wavelet 
vectors  are  the  columns  of  W.  ) 

2.  The  transmitted  waveforms  are  "stretched"  to  twice 
the  W-row-length ,  n,  to  allow  for  2-way  travel. 

An  eight  dimensional  example  with  the  Haar  wavelet  follows  in 
Figure  3.17.  The  example  easily  can  be  extended  to  a  32  or 
64  dimensional  version.  The  Haar  wavelet  is  easy  to 
describe,  but  due  to  it’s  non-smooth  nature  it  is  probably 
not  a  good  wavelet  waveform  for  radar  transmission.  In 
Figure  17,  instead  of  stepping  frequencies  from  one  waveform 
to  the  next,  the  shape  of  the  waveform  changes  in  a  wavelet 
related  fashion.  There  are  8  waveforms  used  to  break  an  8 
foot  long  range  bin  into  8  one  foot  segments  by  wavelet 
transform  (instead  of  FFT).  The  time  duration  of  the 
waveforms  is  16  feet  divided  by  the  speed  of  light. 

There  are  two  speed  advantages  to  transmitting  wavelet- 
related  signals  and  processing  by  wavelets  to  get  fine  range 
prof i les . 

1.  Wavelet  transforms  are  0(n);  FFTs  are  0(nlog2(n)). 

2.  D2 ,  D4,  and  D6  wavelet  transforms  use  1-cycle 
real*real  multiplies;  FFTs  of  I,Q  data  use  6-cycle 
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comp lex*comp lex  multiplies  (4  multiplies, 
and  1  subtract) 


Figure  3.17:  Haar  Waveforms  for  HRRPs 
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4.0  Supplemental  Results 


Figures  4.1  through  4.3  contain  plots  of  partial 
reconstruction  by  D6  wavelets,  of  a  normalized  odd  bounce 
polarization  HRRP  from  a  radar  look  at  the  calibrated 
trihedral  corner  reflector.  Figures  4.4  through  4.6  contain 
similar  plots  for  an  even  bounce  polarization  HRRP  from  a 
radar  look  at  two  dihedrals  that  are  separated  by  about  10 
feet  in  range.  Each  plot  lists  the  number  of  terms  used  in 
the  reconstruction  and  the  normalized  mean  square  error. 

Figures  4.7  through  4.10  show  estimated  target  length 
versus  aspect  angle  for  the  M35  truck.  The  polarization  was 
even  (transmit  left  circular  and  receive  left  circular).  The 
four  plots  are  for  FFT  (Martin  Marietta  method)  processing, 
wavelet  reconstruction  dimension  reduction  by  a  factor  of  2, 
then  4,  then  S  respectively.  Similar  results  are  given  for 
an  odd  polarization  (transmit  right  circular  and  receive  left 
circular)  in  Figures  4.19  through  4.22.  Since  target  length 
will  vary  with  target  aspect,  it  was  decided  that  a  good 
measure  of  length  estimation  stability  would  be  the  standard 
deviation  (over  blocks  of  50  radar  looks)  of  the  length 
estimate.  That  is,  a  small  standard  deviation  would  indicate 
stable  length  estimation.  These  results  are  presented  in 
Figures  4.11  through  4.18  for  an  even  polarization  M35  truck 
and  in  Figures  4.23  through  4.30  for  an  odd  polarization  M35 
truck.  Figures  4.11,  4.13,  4.15,  and  4.17  contain  plots  of 
the  mean  of  the  length  estimates  from  blocks  of  50  M35  looks. 
Figures  4.12,  4.14,  4.16,  and  4.18  contain  plots  of  the 
standard  deviation  of  the  length  estimates  from  blocks  of  50 
M35  looks,  along  with  a  numerical  mean  and  standard  deviation 
of  these  standard  deviations.  This  standard  deviation 
comparison  was  the  figure  of  merit  used  to  determine  the 
stability  of  the  wavelet  methods. 

Figures  4.31  through  4.34  show  estimated  target  length 
versus  aspect  angle  for  the  M60  tank.  The  polarization  was 
even  (transmit  left  circular  and  receive  left  circular).  The 
four  plots  are  for  FFT  (Martin  Marietta  method)  processing, 
wavelet  reconstruction  dimension  reduction  by  a  factor  of  2, 
then  4,  then  8  respectively.  Similar  results  are  given  for 
an  odd  polarization  (transmit  right  circular  and  receive  left 
circular)  in  Figures  4.43  through  4.46.  Target  length 
standard  deviation  versus  aspect  angle  results  for  an  M60 
tank  are  presented  in  Figures  4.35  through  4.42  for  an  even 
polarization  M60  tank  and  in  Figures  4.47  through  4.54  for  an 
odd  polarization  M60  tank.  During  the  analysis  it  was 
discovered  that  a  portion  of  the  data  for  the  M60  tank  was 
found  to  be  unusable.  Although  it  was  too  late  in  the 
analysis  to  determine  the  problem  and  supply  the  correct  data 
to  our  subcontractor,  it  had  no  significant  outcome  on  the 
results . 

For  both  trucks  and  tanks  these  results  show  wavelet 
dimension  reduction  by  a  factor  of  2  or  4  looks  acceptable 
based  on  a  desire  for  small  standard  deviations,  but  when  the 
dimension  is  reduced  by  a  factor  of  8  the  standard  deviation 
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becomes  large  and  erratic.  Consequently,  the  results  show 
that  dimension  can  be  reduced  by  up  to  a  factor  of  4  but  not 
8,  with  no  appreciable  loss  of  stability  in  length  esti¬ 
mation  . 
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Figure  4.2;  Partial  Wavelet  Reconttruction — TVihedral,  Ns(16) 
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Figure  4.7; 


Target  Extent  versus  Aspect  Angle 
Target:  M35,  Threshold:  25,  Polarization:  LL 
Analysis  of  Real  data  (Miutin  Method) 


aspect  angle 

Figure  4.8: 

Target  Extent  versus  Aspect  Angle 
Target;  M35.  Threshold:  25,  Polarization:  LL 
Analysis  of  Real  data  (D6  Level  I) 
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Figure  4.9: 


Target  Extent  versus  Aspect  Angle 
Target  M35,  Threshold:  25,  Polarization:  LL 
Analysis  of  Real  dau  (D6  Level  2) 


aspect  angle 

Figure  4.10: 

Target  Extent  versus  Aspect  Angle 
Target  M35.  Threshold:  25,  Polarization:  LL 
Analysis  of  Real  data  (D6  Level  3) 

aspect  angle 


s.d.  extents  mean  extents 


Figure  4.11 


Average  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.  Target:  M3S.  Threshold:  2S,  Polarization:  LL 
Analysis  of  Real  data  (Martin  Method) 


Block  4  50 


Figure  4.12:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.  Target:  M3S.  Threshold:  2S.  Polarization:  LL 
Analysis  of  Real  data  (Martin  Method) 
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Block  #  50 


mean  =  2.563  s.d.  =  .7725 
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s.d.  extents  level  2  mean  extents  level  2 


Figure  4.15:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.Targeu  M35.  Threshold:  25,  Polarization:  LL 


1  gu  re  4.16:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO,  Target:  M3S,  Threshold  25,  Polarization:  UL 
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mean  extents  level  3 


Figure  4.17:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.  Target*  M3S,  Threshold  25.  Polarization:  LL 


Figure  4.18:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO,  Target*  MSS.  Threshold:  2S.  Polarization:  LL 
2Q^ _ _ _ A^ysis  of  data  (D6  Level  3)  _ 


vehicle  extent  level  1  vehicle  extent 


Figure  4.19:  Target  Extent  versus  Aspect  Angle 
Target  M3S.  Threshold:  23.  Polarization:  LR 
Analysis  of  Real  data  (Martin  Method) 
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Target  Extent  versus  Aspect  Angle 
Target:  M3S,  Threshold:  25,  Polarization:  LR 
Analysis  of  Real  data  (D6  Level  2) 


aspect  angle 


Target  Extent  versus  Aspect  Angle 
Target:  M3S,  Threshold:  25.  Polarization:  LR 
Analysis  of  Real  data  (1)6  Level  3) 
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Figure  4.23:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.Targec  M3S.  Threshold:  25.  Polarization;  LR 


Figure  4.24: 
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Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.Targec  M3S. Threshold:  25. Polarization:  LR 
Analysis  of  Real  data  (Martin  Method) 


s.d.  extents  level  1  mean  extents  level  1 


Figure  4.25:  Average  Target  Extent  versus  Aspect  Angle 

BkxkSize:  SO.  Target;  M3S.  Threshold:  25.  Polarization:  LR 


i  gu  re  4.26:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.  Target:  M3S.  Threshold:  2S.  Polarization:  LR 
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s.d.  extents  level  2  -n  mean  extents  level  2 


Figure  4.27:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO,  Target:  M3S.  Threshold:  2S.  Polarization:  LR 


igure  4.28:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.Targec  M3S.  Threshold:  25,  Polarization:  LR 


mean  =  2.578  s.d,  =  1 .023 
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mean  extents  level  3 


Figure  4.29:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  50.  Target'  M35,  Threshold:  25,  Polarization:  LR 
Analysis  of  Rci  data  (D6  Level  3) 


Block  #  50 


Figure  4.30:  Standard  Deviation  Target  Extent  versus  Asp^t  Angle 
Block  Size:  50,  Target  M35,  Threshold:  25,  Polarization;  LR 
Analysis  of  Re^  data  (D6  Level  3) 
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Figure  4.31 


Target  Extent  versus  Aspect  Angle 
Tnec  M60,T1vtshoId;  2S,Marizxiion:  LL 


vdiicle  extent  level  2 


Figure  4.33:  Target  Extent  versus  Aspect  Angle 
Tagec  M60.11ireshold:  2S.  Ftalvizuioa;  LL 
Anilysii  of  Reel  dsia  (D6  Level  2) 
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Figure  4.34;  Target  Extent  vtrsus  Aspect  Angle 
Target  M60.  Threshold:  2S,n>Iaruaiioa'  LL 
Analysis  of  Real  (laa(D6  Level  3) 
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Figure  4.37:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  50,TirgeL  MdO.llutshold:  2S.Poiaization:  LX. 


F  i  gur  e  4 . 38 ;  Standard  Deviation  Target  Extent  versus  Aspect  Angle 


mean  >  2.955  rd.  s  1.14 
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Figure  4.39:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.  Target:  M60.11ueshol±  2S.Polarizatioa:  LL 
Analysis  of  Real  data  (D6  Level  2) 
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Block  #  50 


Figure  4.40:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.  Target:  M60.  Threshold:  25.Polarizatioa:  LL 
Analysis  of  Real  data  (D6  Levd  2) 
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means 3.417  s.d. >1.114 
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Figure  4.43:  Target  Extent  versus  Aspect  Angle 
Tvgec  M60.ThfeslKdd:  2S.Polarizatiaii:  LR 
Analysis  of  Real  data  (Martin  Method) 


Figure  4.44; 
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Target  Extent  versus  Aspect  Angle 
Target:  M60.  Threshold:  25.Pt>Iarixaiioo:  LR 
Analysis  of  Real  data  (D6  Level  1) 
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Figure  4.45:  Target  Extent  versus  Aspect  Angle 
TvscL*  M60.1hrtsbold:  25,FDliruation:  LR 
Analysis  of  Real  data  ^  Level  2) 
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Figure  4.46:  Target  Extent  versus  Aspect  Angle 
Target  M60,  Threshold  2S.  Polarizatioo:  LR 
Analysis  of  Real  data  ^  Level  3) 
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Figure  4,47:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.Ta^ec  M60.  T1ueshoi±  2S.Mviaiioa:  LR 
_  Aittlyu  of  Real  dm  (Matin  Molnd) 
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Figure  4.48:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  SO.  Target'  M60,  Threshold:  25.  PMauation:  LR 
^ _  Analysis  oTRealdaU  (Martin  Method) 


Figure  4.49:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.Tvgec  M60.Thfeshold:  25.  Poizriiaiion-  LR 
^  Anaiysa  of  Reel  dua(D6  Level  1) 
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Figure  4.50:  Standard  Deviation  Target  Extent  versus  Aspect  Angle 
Block  Size:  S0.Tirgec  M60,  Threshcdd:  TS.PoIaizaiioo:  LR 
_  Anilysis  of  Real  (tea  (D6  Level  1) 
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Figure  4.51:  Average  Target  Extent  versus  Aspect  Angle 

BlockSize:  SO.Tvgee  M60. Uphold:  2S.h>UrixuioD:  LR 
„  Analysis  of  Reel  data  Q>6Levd  2) 
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Figure  4.52:  Standard  Deviation  Tai^et  Extent  versus  Aspect  Ang^e 
BlockSize:  SO.Targec  M60.*nueshold:  2S, Roiaruaiioa'  IR 
2Q _ Analysis  of  Real  data  (D6  Level  2) 
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Figure  4.53:  Average  Target  Extent  versus  Aspect  Angle 

Block  Size:  SO.Tvgec  M60.Thfesh(^‘  2S,  Mvizuion:  LR 
„ _  Analysis  ofReal  (tea  (D6  Level  3) 
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Figure  4.54:  Standard  Deviation  Target  Extent  versus  Aspect  Angie 
Block  Size:  SO.  Target  M60.  Threshold:  2S.h>larizaiioa:  LR 
_  Analysis  of  Real  daa(D6  Level  3) 
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Appendix  A 
List  of  Symbols 


D  :=  Ring  of  dyadic  rational  numbers 
:=  {x  :  X  =  m/2"  ,  m,n€Z,n>0} 

R  :=  Field  of  real  numbers. 


Z  :s  Ring  of  rational  integers. 


W 


Multiplier  for  wavelet  system.  6  Z  and  /i  >  2. 
Wavelet  Coeffident  Matrix  (“WCM^). 


V«o"‘  a?  •••  ON-t  J 


where  N  is  an  integer  multiple  of  the  multiplier  fi. 
If  H  — 2  then  we  write: 


[a] 


Wavdet  Coefficient  Matrix  (“WCM"). 

(oq  ai  •••  \ 

bo  hi  ...  bff^i  j 


N  :=  Number  of  columns  in  a  WCM.  W  is  an  integer  multiple  of  /i. 
If  K  2  then  N  is  even. 
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Appendix  B 

Wavelet  Transform  Theory 


B.l  Wavelets  and  Wavelet  Transforms 

Most  signals  in  science  and  engineering  are  modeled  as  mathematical  func¬ 
tions  for  purposes  of  analysis.  In  order  to  separate  or  examine  certain  impor¬ 
tant  features  or  characteristics  of  the  signal,  the  function  is  often  expanded 
in  terms  of  basis  functions  that  span  the  space  or  a  subspace  that  the  sig¬ 
nals  of  interest  reside  in.  The  most  common  example  of  this  is  the  Fourier 
transform  where  a  signal  that  oripnates  in  the  time  domain  is  reformulated 
in  the  frequency  domain  by  expanding  the  function  in  terms  of  trigonometric 
or  complex  exponential  basis  functions.  This  basis  is  most  appropriate  when 
the  signals  have  periodic  components  or  are  produced  by  systems  that  are 
modeled  by  constant  coefficient  differential  or  difference  equations. 

Consider  a  periodic,  possibly  complex-valued,  signal  g{t)  that  is  square 
integrable  over  the  range  (0  <  t  <  1}  and  with  period  one  so  that 

9{t)  =  9{t  +  1).  (B.l) 

This  function  can  be  expanded  in  a  Fourier  series  of  the  form 

,(!)=  f;  (B.2) 

with  the  coefficients  given  by 

4.= /',(()  (B.3) 
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which  is  an  inner  product  of  g{t)  with  the  basis  functions.  Similarly,  one 
defines  the  wavelet  transform  with  respect  to  a  basis  of  wavelet  functions. 

The  wavelet  basis  is  generated  by  h  ft  x  N  matrix  [a]  ,  where  /i  and  N 
are  positive  integers  and  N  ia  k  multiple  of  /i .  The  matrix 


[a]  := 


"o. 


V  oo"* 


“atm  / 


(B.4) 


is  a  wavelet  coefficient  matrix  (‘‘WCM”)  if  it  satisfies  the  scaling  conditions 


N-t 


(B.5) 


E4  =  /'<o.-  (B.6) 

k 

where  equals  1  if  t  =  j  and  0  otherwise.  The  overbar  denotes  complex 
conjugation  and  /  is  an  integer.  The  sums  over  k  axe  finite  sums  since  only 
finitely  many  of  the  numbers  al  are  different  from  zero. 

The  positive  integer  ft  is  called  the  multiplier  of  the  wavelet  system  and 
N  is  called  its  length; 

This  matrix  of  numbers  provides  coefficients  for  the  vector  of  recursions 
V7'[a](x)  =  ^  aiv?°[a](2x  -  k)  (B.7) 

ibsO 

which  implicitly  define  the  wavelet  scaling  function  and  explicitly  define 
the  basic  wavelet  functions  ^'[a](x),  1  <  i  <  p-  Observe  that  only  ^°[a] 

appears  on  the  right  hand  side.  The  functions  ^*[a](x)  are  defined  for  all  real 
numbers  x  €  R< 

The  fundamental  fact  about  systems  of  compactly  supported  wavelets  is 
that  the  collection  of  functions 

Basis[a]  :=  —  it) : 

0<»</i,  Z} 
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form  a  basis  for  L*  spaces,  in  particular  for  X’(R)- 

We  now  focus  on  the  case  of  ft  =  2  and  will  use  the  following  simplified 
notation  for  the  scaling  function  ^(O  =  9°[a](t)  and  the  wavelet  function 
=  ^*[a](t)  where  t  €  R*  It  can  be  shown  that  if  the  coefficients  of  this 
equation  satisfy  the  wavelet  conditions,  stated  above,  the  solution  ^(t)  will 
be  orthogonal  to  integer  translates  of  itself  and  can  be  normalized  such  that 

<  V»(0.  (fi{t-k)>  =  j  v>(0  -k)di-  Sojc  (B.8) 

This  means  the  set  of  basis  functions 

w(l)  =  vKf  - ')  (B.9) 

spans  a  subspace  Vo  in  and  the  coefficients  of  an  expansion  within  this 
subspace  can  be  calculated  as  simple  inner  products.  The  feature  of  scaling 
functions  that  makes  them  attractive  for  signal  processing  is  their  ability  to 
model  signal  properties  that  are  related  to  the  independent  variable  t.  One 
can  increase  the  size  of  the  subspace  spanned  by  the  scaling  functions  by 
using  —  k)  which  spans  a  subspace  V,-.  One  can  show 

that  Vo  C  Vj  C  Vj  C  '  • 

The  features  of  a  signal  can  often  be  better  described  by  defining  a  slightly 
different  set  of  orthogonal  basis  functions  that  span  the  differences  between 
the  spaces  spanned  by  the  various  scales  of  the  scaling  function.  These  new 
functions  are  the  wavelets.  The  basic  wavelet  is  defined  in  terms  of  the  scaling 
function  by 

m  =  -  k) .  (B.IO) 

* 

It  b  the  prototype  of  a  class  of  orthonormal  basis  functions  of  the  form 

-  it)  (B.ll) 

where  2'  b  the  scaling  of  t.  2~^k  is  the  translation  in  f,  and  2^^^  maintains 
the  unity  norm  of  the  wavelet.  We  shall  say  that  j  b  the  base-2  logarithm 
of  the  s^e.  If  Wj  is  the  subspace  of  L*  —  ^’(R)  spanned  by  the  integer 
translates  of  the  wavelet  2^^’0(2^t),  additional  disjoint  subspaces  are  spanned 
by  integer  translates  of  the  wavelets  for  each  different  scale  index  j,  that  b, 
by  the  functions  ^jjt{t).  The  relationship  of  the  various  subspaces  can  be 
seen  from  the  following  expressions: 

VoC  Vi  CVaC  =  /i’(R),  (B.12) 
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Vo©Wo  =  V,  ,  (B.13) 

V,_i  ©  W^-,  =  V>  ,  (B.14) 

L2  =  Vo©Wo©Wi©---  (B.15) 

and  indeed,  if  we  allow  j  to  run  over  all  integers,  then 

I*  =  . . .  W_2  ©  w_,  ©  Wo  ©  W, . . .  Wj . . .  .  (B.16) 

This  states  that  the  set  of  basis  functions  formed  from  and  span 

all  of  L*  and,  therefore,  any  function  in  can  be  written 

<7(0=  £  c<V’'(0  +  £  £  (B.17) 

/=— oo  j=0k=—oo 

with  the  coefficients  expressed  by 

Cl  =  j  9{t)  <pi{t)  dt  (B.18) 

and 

dj,k  =  [  9{t)  ^j,k{t)  dt.  (B.19) 


The  basis  functions  9>i{t)  and  0jfc(O  are  numerical  valued  functions  of  nu¬ 
merical  variables;  they  have  no  physical  dimension.  In  the  expansion  formula 
(B.17)  the  argument  2H  —  k  oi  rjf  \s  a  pure  number  so,  writing 

the  quantity  2^  has  the  dimension  i.e.  frequency,  and  fc/2^  has  the 
dimension  time. 

These  wavelet  coefficients  completely  and  uniquely  describe  the  original 
signal  and  can  be  used  to  represent  it  in  a  way  similar  to  Fourier  coefficients. 
Because  of  the  orthonormality  of  the  basis  functions,  there  is  a  version  of 
Parseval’s  theorem  that  relates  the  energy  of  the  signal  9{t)  to  the  energy  in 
each  of  the  wavelet  expansion  components  and  their  wavelet  coefficients  by 
the  formula  ^ 

llj||’  =  EW'  +  E  E 

I  j>0k=—oo 
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This  is  one  reason  why  the  orthonormality  is  so  important.  Daubechies  [2] 
showed  that  the  translates  of  the  scaling  function  and  the  translated  dilations 
of  the  wavelets  are  orthonormal,  and  all  of  the  these  functions  have  compact 
support  (i.e.  are  non-zero  only  over  a  finite  region)  if  there  are  only  a  finite 
number  of  non-zero  coefficients  C),  in  the  recursive  scaling  equation  (B.7). 
This  provides  the  time  localization  that  is  particularly  desirable  for  analyzing 
both  the  time  and  the  frequency  behavior  of  transient  signals. 

Note  that  there  is  an  infinite  set  of  scaling  functions  and  wavelets  that 
can  be  obtained  by  choosing  different  coefficients  an  in  (B.7). 


B.2  Wavelet  Approximation  to  the  Fourier 
Transform 

In  this  chapter  we  state  the  relationships  between  the  wavelet  transform  and 
the  Fourier  transform  as  derived  by  Resnikoff  and  Bumis  [5]. 

In  this  section  we  restrict  attention  to  periodic  functions  g{t).  Let  mbox^^(t+ 
1)  =  g(^t).  Then  the  scaling  function  and  wavelet  coefficients  in  (B.18)  and 
(B.19)  are  also  periodic.  The  scaling  function  coefficients  c(I)  calculated  from 
(B.18)  and  using  (B.l)  are 

cj  =  /  (^(i-O  g(t)  dt=  f  (p{t)  g{T  +  l)dT=  I  g^r)  dr  =  Cq. 

(B.21) 

Using  the  ‘‘partition  of  unity”  property  of  scaling  functions  which  states 

f;v(f-0  =  l  {B.22) 

im—oo 

we  can  write  (B.21)  as 

^  =  53  j[  ¥»(’■-  0  ^(i")  =  j[  53  V(»‘  “  0  ‘i’’  =  dr  =  c. 

(B.23) 

Therefore,  C(  is  a  constant  and  not  a  function  of  /.  Note  the  interesting 
property  of  (p(t)  which  causes  the  scaling  function  coefficient  of  a  periodic 
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function  to  be  equal  to  the  DC  Fourier  series  coefficient  for  any  scaling  func¬ 
tion  satisfying  (B.7); 

r  V>(0  y(0  *  =  r  ^(0  dt .  (B.24) 

J—oo  Jo 

Using  the  partition  of  unity  property  (B.22)  again  shows  the  scaling  function 
term  of  the  expansion  in  (B.17)  is  a  constant. 

=  =  (B.25) 

/  I 

In  a  like  manner,  using  (B.19)  we  calculate  the  wavelet  coefficients  in  the 
expansion  (B.17)  by 

<<,.»+*•  =  2''" / HVl  -k-K) g(t)  di  =  V/’  j - 2-’K)  - k)  j(l)  * 

(B.26) 

which,  after  the  change  of  variables  r  =  <  —  2^K  and  /  =  2~*K  or  /if  =  2^7, 
diMK  =  V/’  J  0(2^t  -  k)  ,(r  +  2-iK)dT  =  djj..  (B.27) 

This  states 

djjt  =  (B.28) 

or,  the  wavelet  coefficients  dy,*  are  periodic  in  k  with  period  2^.  We  now 

use  these  properties  with  the  basic  definition  of  the  Fourier  series  coefficients 
from  (B.3). 

6„=  (B.29) 

Jq 

Substituting  the  wavelet  expansion  from  (B.17)  for  g{t)  into  (B.29)  gives 

f  c  ifn  =  0; 

K^\  (B.30) 

I  ra"5«0 

where  0(/)  is  the  Fourier  transform  of  0(t)  defined  by 

^-00 
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with  an  inverse  given  by 


v.(o  =  r  «/) 

«/— OP 


(B.32) 


An  alternative  form,  which  uses  the  Fourier  transform  of  the  basic  wavelet 
0(t)  rather  than  0j>(Oi  [^]* 


if  n  =  0; 


ifnjiO 


(B.33) 


where  dj{f)  is  the  discrete  Fourier  transform  (DFT)  of  for  various  scales 
j  and  is  defined  by  [1] 


k=0 


(B.34) 


with  an  inverse  given  by 


i7wtnk/N 


(B.35) 


The  formula  for  expressing  the  wavelet  coefficients  in  terms  of  the  Fourier 
coefficients  is  easily  obtained.  From  (B.19)  we  have 


djjk  =  /  y(0  0i.*(O  dt. 

OP 


which,  after  using  (B.2)  and  reordering,  gives 


<*»>  =  E  Jw  r  * 

NS— CO  OO 


(B.36) 


Yd 

ns— OP 

and,  from  (B.24),  the  scaling  function  coefficient  is  given  by 

c*  =  ho  =  c. 


(B.37) 


(B.38) 
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An  alternative  derivation  of  (B.37)  is 

(B.39) 

where  the  discrete-time  Fourier  transform  is  defined  by 

DTFT{i.}  =  i(/)  =  f; 

(B.40) 

and  the  inverse  is 

i.=  / 

Jt\ 

(B.41) 

These  are  the  basic  relationships  between  the  Fourier  transform  and  wave¬ 
let  transform.  Equation  (B.33)  relates  the  Fourier  series  expansion  coeffi¬ 
cients  of  ^(t)  to  the  Fourier  transform  of  the  fundamental  wavelet  and  the 
DFT  of  the  wavelet  expansion  coefficients  of  g{t).  Equation  (B.39)  expresses 
the  wavelet  expansion  coefficients  of  g{t)  in  terms  of  samples  of  the  DTFT 
of  the  product  of  samples  of  the  Fourier  transform  of  the  basic  wavelet  and 
the  Fourier  series  coefficients. 

In  practice,  a  band  limited  signal  sampled  at  the  appropriate  rate  implies 
a  finite  sum  in  (B.33)  that  corresponds  to  the  finest  scale  of  j  sz  J  and  the 
Fourier  series  is  likewise  finite.  The  independent  variable  /  is  frequency  in 
Hertz  or  cycles  per  second,  d(m)  is  periodic  with  period  N  =  2\  and  x{f)  is 
periodic  with  period  one.  Clearly,  the  DFT  or  its  inverse  can  be  calculated 
with  the  FFT  [1]. 

Energy  and  Parseval*s  Theorem 

Orthonormal  basis  systems  allow  direct  calculation  and  interpretation 
the  energy  in  a  signal  partitioned  in  both  the  time  and  the  expansion 
domains.  Parseval’s  theorem  for  the  Fourier  series  (B.2)  states 

E  w’.  (B.42) 

fiS-*00 

The  ‘‘power”  in  a  signal  is  proportional  to  the  square  of  the  signal  (e.g. 
voltage,  current,  force,  or  velocity)  and,  therefore,  the  energy  is  given  by 
the  integral  of  the  square  of  the  signal  magnitude.  Parseval’s  theorem  states 
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how  the  total  energy  is  partitioned  in  the  frequency  domain  in  terms  of  the 
partition  provided  by  the  the  orthonormal  basis  functions.  For  the  general 
wavelet  expansion  of  (B.17),  Parseval's  theorem  is 

/  IsWl’ <«  =  f:  Icp  +  E  E  (B.43) 

with  the  energy  in  the  expansion  domain  partitioned  in  time  by  /  and  k  and 
in  scale  by  j.  For  the  case  of  periodic  functions,  the  relationship  reduces  to 

('  ls(i)l’  -ft  =  Icl’ + E  E  ■  (B.44) 

One  can  show  that  |^(/)|  — »  0  as  /  — ♦  0  and  also  as  /  -♦  oo.  Therefore, 
there  will  be  a  band  of  frequencies  where  most  of  the  energy  in  ^(/)  is 
concentrated.  Likewise,  for  many  signals,  the  energy  will  be  concentrated  in 
a  region  of  the  (j,  k)  plane.  Because  of  this  concentration  and  using  (B.33), 
the  energy  of  the  signal  ^(t)  at  frequency  n  and  at  scale  j  and  time  k  is 
approximately  measured  by 

(B.45) 

If  most  of  the  energy  in  ^(/)  occurs  around  frequency  /o,  then  /o  =  n/2^ 
relates  the  dominant  Fourier  frequency  /o  to  the  dominant  wavelet  scale  j. 
Scale  and  frequency  are  independent  primitive  concepts,  but  the  selection 
of  a  wavelet  basis  establishes  a  connection  between  them,  and  the  results  of 
this  chapter  allow  one  to  move  between  the  two  descriptions,  using  the  one 
most  appropriate  for  a  particular  problem.  The  simple  partitioning  of  the 
energy  content  of  a  signal  among  frequencies  has  been  generalized  to  include 
the  parameters  of  time  and  scale.  In  addition,  we  have  at  our  disposal  the 
choice  of  wavelet  systems,  which  is  controlled  by  the  choice  of  ai,  in  (B.7), 
and  determines  the  detailed  nature  of  the  relationship  between  frequency  and 
scale. 
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B.3  Computational  Complexity  of  the  Wave¬ 
let  Approximation  to  the  Fourier  Trans¬ 
form 

This  section  is  concerned  with  the  computational  efficiency  of  wavelet-based 
approximations  to  the  Fourier  IVansform.  The  computational  cost  for  several 
wavelet  approximations  of  the  Fourier  transform  are  presented  and  compar¬ 
isons  are  made  with  the  computational  complexity  of  the  FFT. 

In  the  most  applications,  the  size  of  the  data  set  is  the  determining  factor. 
We  will  refer  to  a  technique's  asymptotic  complexity  as  representative  of 
how  quickly  the  number  of  computations  grows  as  the  problem  size,  AT,  in¬ 
creases.  Comparisons  with  the  FFT  are  natural  here,  since  the  FFT  is  a 
commonly  used  and  well  understood  algorithm.  We  consider  wavelet  trans¬ 
forms  which  resolve  to  sub-band  frequency  components  in  part  or  perhaps 
all  of  the  signal  bandwidth.  No  attempt  will  be  made  to  assign  a  value  to  a 
particular  scheme,  since  the  level  of  approximation  appropriate  for  a  given 
application  can  only  be  determined  by  extensive  experiments.  The  quality  of 
the  approximation  varies  from  approach  to  approach,  and  many  simple  and 
consistent  optimizations  can  be  made,  but  usually  require  specialization  to 
a  particular  set  of  circumstances. 

The  simplest,  and  also  least  accurate  technique  for  estimating  the  DFT 
of  a  signal  is  a  direct,  balanced  tree  approach  with  filters  optimized  for  sharp 
out-of-band  rejection.  Typically  the  Daubechies  wavelets  are  used  for  this 
approach,  although  slight  improvements  in  the  quality  of  the  results  are  pos¬ 
sible  with  other  choices.*  The  computational  complexity  of  pure  wavelet 
methods  are  analyzed  in  section  2.5. 

Additional  improvements  in  the  quality  of  the  results  and  the  speed  of 
computation  is  possible  when  full  frequency  resolution  within  a  limited  sub- 

*Ileal  wavelet  filters  have  spectral  characteristics  which  are  coqjugate-synunetric,  aad 
these  would  combine  large  and  small  frequencies  in  the  same  band.  To  avoid  this  aliasing, 
a  modification  was  made  to  the  analysis  and  synthesis  filters.  Simply  put,  the  filter 
characteristics  were  "rotated”  so  that,  for  a  full  bandwidth  complex  signal,  whose  Fourier 
transform  has  the  full  2N  degrees  of  freedom,  the  low-pass  side  of  a  Mgle  analysis  stage 
captures  the  lower  half  of  the  spectrum  and  the  high  pass  the  upper.  In  the  real  case, 
the  Fourier  transform  is  conjugate-symmetric,  the  "full”  bandwidth  anuunts  to  only  ft/2, 
half  the  sample  rate,  and  "low-frequency”  is  then  interpreted  as  being  /«/4. 
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band  is  the  desired  result.  This  is  accomplished  by  using  the  wavelet  trans¬ 
form  to  separate  sub-bands,  and  then  using  an  FFT  on  the  data  within  the 
band  to  reconstruct  a  ‘‘local’*  approximation  to  the  “full”  Fourier  transform. 
This  can  be  done  either  (1)  solely  within  that  band,  or  (2)  throughout  the 
entire  bandwidth.  If  the  sub-bands  used  in  such  a  reconstruction  capture 
a  significant  fraction  of  the  total  energy  contained  in  the  signal,  excellent 
approximations  can  be  constructed  quite  cheaply.  Notably,  the  locations  of 
peaks  within  bands  are  properly  located  and  reconstruction  away  from  edges 
of  sub-bands  is  quite  accurate.  The  complexity  of  such  a  technique  is: 

OPS  =  a{K)N(l  -R)  +  NR(3  +  log(iV/?))  +  RN  (B.46) 

where 

•  Af  =  the  number  of  data  points*, 

•  K  ^  the  size  of  the  wavelet  filter  (assumed  to  be  constant  throughout 
the  procedure; 

•  R  =  the  portion  of  the  bandwidth  of  the  signal  which  is  resolved  (as¬ 
sumed  to  be  1/2^  for  some  L). 

For  this  estimate  we  assume  that  the  reconstruction  is  done  purely  in-band, 
and  only  a  single  band  is  used  (i.e.,  we  don’t  expand  further  than  necessary 
to  locate  the  energy,  since  that  is  wasted  work  if  we  truly  want  the  DFT). 

Figure  B.2  shows  the  results  of  this  calculation  superimposed  on  the  re¬ 
sults  for  a  straight  wavelet  approach.  The  tightly  clustered  lines  represent 
the  hybrid  case,  and  clearly  show  a  usable  advantage  (i.e.,  reasonable  length 
filters  are  justified)  over  both  the  FFT  and  the  straight  wavelet  methods  for 
problem  sizes  starting  around  64  data  points. 

In  addition  to  the  results  stated  there  we  note  that  significant  improve¬ 
ments  may  be  made  by  combining  these  methods  with  phase-unwrapping  in 
order  to  move  the  “mass”  of  the  return  to  the  region  of  best  filter  perfor¬ 
mance,  namely  zero  frequency  for  real  filters  and  /  =  /•/4  for  the  “rotated” 
filters. 

The  approach  of  relocating  the  center  of  the  “target”  frequency  band  is 
the  essence  of  the  phase  unwrapping  technique,  where  the  centroid  is  moved 
to  the  center  of  the  DC  or  zero  frequency  band.  Real  Daubechies  wavelets 
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Figure  B.2:  Computational  Complexity  of  the  Wavelet  Transform  vs.  the 
FFT:  Filter  Length  vs.  N  for  equivalent  complexity  for  various  fractions  of 
the  total  bandwidth  (R)  with  the  addition  of  the  case  when  the  FFT  is  used 
to  perform  the  final  frequency  analysis  in  a  subband 

have  optimal  performance  at  D.C  and  minimum  aliasing.  Phase  unwrap¬ 
ping  carries  an  overhead  of  N  complex  multiplies  per  N  input  points.  In 
addition,  calculation  of  the  FFT  for  in-band  components  takes  R  additional 
operations  while  calculation  of  the  entire  FFT  from  inband  components  takes 
appronmately  N  additional  computations.  Ndther  of  these  has  a  significant 
contribution  to  the  total  computation  for  large  problem  sizes,  since  both  are 
overwhelmed  by  the  O(nlog3(n))  behavior  of  the  DFT’s.  The  number  of 
operatbns  is  then: 

OPS  ^N  +  RN-\- aN{l  -R)  +  NR{3  +  log(AfR))  (B.47) 
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Figure  B.3:  Computational  Complexity  of  the  Wavelet  Transform  with  Phase 
Unwrapping  vs.  the  FFT:  Filter  Length  vs.  N  for  equivalent  complexity  for 
various  values  of  R  for  the  case  of  using  phase  unwrapping  and  the  FFT  to 
perform  the  final  frequency  analysis  in  a  subband 


Figure  B.3  shows  the  trade-off  between  this  method  and  an  FFT  on  the 
full  problem,  reflected  in  the  allowable  filter  size  for  equivalent  quantities  of 
computations. 
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Appendix  C 

Phase  Unwrapping 


In  this  section  we  describe  a  method  for  preprocessing  MMW  radar  data  to 
concentrate  the  target  energy  into  the  largest  scale  wavelet  coefficients.  It 
is  called  phase  unwrapping.  It  can  be  used  to  process  radar  fine  range  data 
from  phase  history  data  sequences  whose  bandwidth  is  less  than  one  half  the 
sampling  rate. 

This  section  has  the  following  objectives: 

•  to  describe  the  phase  unwrapping  algorithm 

•  to  describe  a  method  of  radar  imaging  based  on  combining  phase  un¬ 
wrapping  with  the  one-sided  *‘Maliat”  algorithm 

•  to  present  numerical  examples  using  actual  radar  phase  history  data 

The  phase  unwrapping  algorithm  is  based  on  the  topological  concept  of 
winding  number.  Let  z  :  [a,6]  — »  C  be  a  differentiable  complex  valued 
function.  Define  the  winding  number  of  z  over  the  interval  [a,  6]  by 

W(z^a^b)si  f  Real{w{t)/z{t))dtyW  dz/idt.  (C.l) 

Jtmm 

Hence  H^(z,a,6)  represents  the  amount  of  angular  winding  about  0  (in 
the  counterclockwise  direction)  of  the  point  z(t)  as  t  moves  from  a  to  b.  If  z(t) 
represents  (interpolated)  phase  history  data  formed  from  a  rough  reflector 
(several  scatterers  per  resolution  cell)  then  it  is  generally  understood  that 
z(t)  may  be  approximately  modelled  by  a  complex  Gaussian  random  process 
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Figure  C.l:  Radar  Centroids  for  M60  Tank  Target 

having  power  spectral  density  function  I{y).  Furthermore,  I  represents  the 
expected  value  of  |z(y)P  where  z  denotes  the  Fourier  transform  of  z. 

The  expected  value  E(W)  of  W(z,  a,  b)  equals  2x(5— a)  times  the  centroid 
of  /,  thus 

E{W)  =  2x(6  -a)J  /(x)x  dx/  J  /(x)  dx.  (C.2) 

The  factor  2t  arises  because  of  its  appearance  in  the  definition  of  the  Fourier 
transform. 

The  radar  centroids  and  phase  centers  were  calculated  for  2000  different 
aspect  angles  of  the  M60  tank  target  to  evaluate  the  stability  of  the  phase  un¬ 
wrapping  algorithm.  Figure  C.l  plots  the  results  of  the  centroid  calculations 
and  figure  C.2  plots  the  results  of  the  phase  unwrapping  method.  Statistical 
analysis  of  these  two  sets  of  data  reveals  that  the  mean  difference  is  0.33 
range  cells  and  the  standard  deviation  of  the  difference  is  1.32  indicating  a 
slight  bias  but  otherwise  very  consistent  results.  The  correlation  coefficient  r 
between  the  two  sets  of  data  is  0.9986  indicating  that  the  phase  unwrapping 
method  is  a  near  perfect  predictor  of  the  true  radar  centroid. 

The  phase  unwrapping  algorithm  consists  of  two  parts: 

1.  estimate  the  total  phase  winding  of  a  phase  history  function  z{t)  from 
its  sample  values  Zk, 

2.  remove  this  phase  winding  by  multiplying  Zk  by  an  appropriate  complex 
sinusoidal  sequence. 
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Figure  C.2:  Phase  Centers  Calculated  by  Unwrapping  Method  for  M60  Tank 
Target 
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Figure  C.3:  Raw  Range  Data  and  DFT  for  Trihedral  Target 


Figure  C.4:  Phase  Unwrapped  Data  and  DFT  for  Trihedral  Target 


90 


Figure  C.5:  Raw  Range  Data  and  DFT  for  M60  Tank  Target 
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Figure  C.6:  Phase  Unwrapped  Data  and  DFT  for  M60  Tank  Target 

The  first  step  requires  estimating  W (s, a,  b)  =  W{z,k,k  +  \)  from 

sample  values  zi,  under  the  assumption  that  the  bandwidth  of  z(x)  is  less 
than  half  than  the  sampling  rate.  In  this  case  either  Zk  or  (— l)*Zik  is  a  low 
frequency  content  signal  dk.  Let  an  denote  the  angle  change  between  dk+i  and 
dk.  Then  the  average  magnitude  of  at,  is  less  than  90  degrees.  This  criteria 
can  be  used  to  select  di,  =  zu  or  dk  =  (— Furthermore,  ak  (nearly 
always)  represents  either  W{z^  k^k+l)  or  — t/2  +IF(z,  it,  ^  +  1).  Therefore, 
VF(z,a,5)  else  W{z,a,b)  w  (6  -  a)T/2  +  a*.  Let 

W  denote  the  total  phase  winding. 

The  second  step  consists  of  multiplying  z*  by  the  sequence  to 

obtain  s*  =  z*e”***^/^*“*l.  By  the  assumptions  on  z*  and  by  (C.2),  s*  has  a 
low  frequency  content. 

To  form  a  radar  image  from  Zk  apply  the  following  steps: 


NonnaUzed  Mean  Squared  Error 


Number  of  Psi  Terms 

Figure  C.7:  Comparison  of  Approximation  Error  for  Trihedral  Target  for 
Original  and  Phase  Unwrapped  Data 


1.  estimate  W  from  Zk, 

2.  compute  si,  =  r*e”***7(*-«)^ 

3.  apply  the  one-sided  recursion  method  to  Sk  to  form  a  radar  image, 

4.  if  desired,  translate  this  image  HY(2r)  pixels  to  the  right  to  locate  the 
image  in  the  correct  position. 

These  three  steps  yield  an  approximation  of  the  radar  image  that  is  formed 
by  computing  the  (Discrete)  Fourier  transform  of  Zk- 

Figures  C.3  and  C.5  are  plots  of  the  (raw)  real  components  and  DFT  for 
a  calibration  target  and  the  M60  tank  target.  Note  the  oscillatory  behavior 
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Normalized  Mean  Squared  Error 


Figure  C.8:  comparison  of  Approximation  Error  for  Dihedral  Target  for  Orig¬ 
inal  and  Phase  Unwrapped  Data 


j 


of  the  raw  radar  data.  The  results  of  "phase  unwrrapping”  the  radar  data  for 
these  targets  are  presented  in  figures  C.4  and  C.6.  In  these  figures  note  the 
dramatic  reduction  of  oscillatory  behavior  which  allows  these  function  to  be 
well  represented  by  only  large  scale  wavelet  terms.  There  is  essentially  no 
loss  of  shape  information  caused  by  the  phase  unwrapping  process  indicating 
that  it  can  be  used  to  reduce  the  computational  complexity  without  adding 
any  significant  artifacts. 
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Appendix  D 

Martin  Marietta  Radar  Data 
and  Pulse  Compression 
Processing 


This  section  examines  the  problem  of  radar  range  profile  processing  that 
may  be  amenable  to  solution  by  employing  wavelet  analysis  in  addition  to  or 
in  place  of  conventional  Fourier*based  analytical  methods.  The  unclassified 
radar  detection  processing  problem  is  discussed  in  detail  by  Einstein  [3]  and 
Skolnik  [6]. 

The  problem  of  obtuning  sufficient  target  signal-to-noise  ratio  for  MMW 
radar  systems  with  low  transmitter  power  is  solved  by  increasing  the  duration 
of  the  transmitted  waveforms,  thereby  increasing  the  transmitted  energy. 
Conventional  ^pulsed"  radar  systems  have  their  range  resolution  determined 
by  the  pulsewidth  of  the  transmitted  waveform  and  increasing  the  pulsewidth 
reduces  the  range  resolution.  This  limitation  is  overcome  by  coding  the 
transmitted  waveform  such  that  processing  (i.e.,  matched  filtering)  in  the 
receiver  can  recover  compressed  pulse  information  with  an  equivalent  pulse 
duration  approximately  equal  to  the  redprocal  of  the  spectral  bandwidth 
of  the  modulated  long  pulse  duration  waveform.  There  are  several  pulse 
compression  methods  available  (see  [6],  §10)  however,  here  we  will  focus  on 
the  problem  of  processing  time-frequencg-coded  waveforms. 
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D.l  Time-Frequency  Pulse  Compression 

A  range  finding  radar  transmitter  emits  frequency  pulses  of  electromagnetic 
energy  of  the  form^ 


where  fk  is  the  (selectable)  constant  frequency,  Ok  is  an  arbitrary  phase  angle, 
and  ^*(0  represents  the  phase  drift  of  the  oscillator  with  respect  to  the 
constant  frequency  /a  as  a  function  of  time.  It  is  conventional  to  select  the 
time  origin  to  coincide  with  the  time  of  emission  of  the  pulse.  The  emitted 
frequency  will  be  assumed  to  serve  as  the  standard  reference  for  the  return 
signal.  By  definition, 

*(0)  =  0  and  J^»(0)  =  0. 

The  signal  reflected  from  a  point  target  is  the  original  transmitted  signal 
delayed  by  the  transit  time  and  diminished  by  spherical  spreading  and  the 
complex  reflection  coefficient  of  the  target  point.  Thus  the  return  signal  has 
the  form 

y{t)  —  i4fu(t  —  a:)e^^‘  (D*2) 

where  the  targets  complex  reflection  coefficient  is 


and  X  is  the  round  trip  time  from  the  radar  to  the  target.  The  range  r  of  the 
target  b 


r  =  cx/2  (D.3) 

where  c  denotes  the  speed  of  light. 

Upon  reception,  the  reflected  signal  is  mixed  with  the  complex  conjugate 
of  the  emitted  signal  to  provide 

«(<)  s  y{t)u{i)  =  -  z)  x  ^(<)  =  (9.4) 

where  *‘u(t)”  denotes  the  complex  conjugate  of  u{t). 

*  V"  denotes  y/^. 
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A  pulse  of  finite  duration  will  attenuate  the  return:  if  the  pulse  shape 
function  hit-*  5(t),  then  the  return  will  be  multiplied  by  S(t  —  x). 

In  accordance  with  Huyghens*  principle,  a  distributed  target  will  produce 
a  return  that  is  the  superposition  of  the  point  returns.  If  p(x)  represents  the 
complex  scattering  profile  (i.e.  coefficient)  of  the  target  and  s(x)  represents 
the  transmitted  pulse  shape,  then  the  received  signal  due  to  a  differential 
scattering  element  p(x)dx  located  at  range  delay  time  x  will  be 

dE{fk,  <)  =  p(x)s(<  -  (D.5) 


and  the  received  voltage  (prior  to  mixing  with  the  reference  signal)  is 


E(fk,  t)  =  r  p(x)s{t  -  .  (D.6) 

J—OO 

Einstein  [3]  shows  that  if>{t  —  x)  ^  0.  With  this  assumption,  and  after 
mixing  with  the  reference  signal,  the  complex  voltage  at  the  output  of  the 
detector  becomes 

V(/*,0  =  r  p(x)s(<  -x)c-"*’'^‘*dx  . 

J^oo 

Define 

Q(M= 

^—00 

Then  the  Fourier  transform 

q{x,  t)  :=  Qifk,  <)  =  p{x)s{t  -  x)  (D.9) 

represents  the  scattering  profile  p(x)  of  the  target,  weighted  by  the  pulse 
shape  s(t  —  x). 

Radar  returns  from  stepped  frequency  pulse  trains  may  be  processed  ei¬ 
ther  coherently  or  noncoherently.  Coherent  processing  consists  of  calculating 
the  Fourier  transform  of  the  N  coherently  detected  complex  pulse  returns 
from  each  range  cell,  where  N  denotes  the  number  of  pulses  in  the  train. 
The  output  is  a  “high  resolution”  range  profile  of  the  scatterers  within  the 
range  cell. 

Noncoherent  processing  consists  of  calculating  the  Fourier  transform  of 
the  square  of  the  magnitude  of  the  of  the  pulse  returns  from  a  given  range 


(D.7) 

(D.8) 
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cell,  and  produces  the  autocorrelation  of  the  range  profile.  In  both  cases 
the  range  resolution  is  detormined  by  the  total  frequency  spread  of  the  pulse 
train. 

Coherent  Processing 

For  coherent  processing  q(x,t)  is  realized  as  the  Fourier  transform  of 

i.e. 

,(!,()=  r  (D.IO) 

oo 

This  is  approximated  by  (  from  D.6) 

p{x)s(t  -  x)  «  /“  Vif,  .  (D.l  1) 

OO 

assuming  that  the  phase  drift  V’ibCO  oscillator  is  small. 

When  using  a  frequency-diverse  pulse  train  the  measurements  of  V{f,i) 
are  not  made  over  a  continuum  of  frequencies  but  only  for  a  discrete  set  of 
frequencies  {fk)-  Thus  the  continuous  transform  must  be  approximated  by 
a  discrete  (and  finite)  transform. 

For  this  case  Einstein  [3]  shows  that 

“The  Discrete  Fourier  TVansform  (DFT)  is  a  matched  filter 
for  the  returns  from  a  frequency-diverse  pulse  train,  under  the 
following  conditions: 

1.  Each  of  the  pulse  returns  is  coherently  detected  by  using  the 
transmitted  carrier  for  a  given  pulse  as  the  coherent  reference 
for  that  pulse.  Thus,  a  different  coherent  reference  is  used 
for  each  pulse  in  the  train. 

2.  The  spadng  between  the  different  frequencies  at  which  the 
pulses  in  the  train  are  transmitted  is  uniform. 

3.  There  is  no  relative  radial  motion  between  the  target  and 
the  radar  during  the  transmission  of  the  pulse  train. 

4.  There  are  no  nonlinear  differential  phase  shifts  through  the 
radar  system  between  the  different  carrier  frequencies.” 
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Item  #2  is  required  for  application  of  the  DFT  to  recover  a  time  domain 
signal  from  the  frequency  return  signal  data.  This  applies  to  both  coherent 
and  noncoherent  processing. 

Noncoherent  Processing 

Noncoherent  processing  consists  of  calculating  the  Fourier  transform  of 
the  square  of  the  magnitude  of  the  of  the  pulse  returns  from  a  given  range 
cell,  and  produces  the  autocorrelation  of  the  range  profile.  Noncoherent 
processing  offers  the  following  advantages  and  disadvantages  compared  with 
coherent  methods: 

•  The  transmitter  frequency  stability  requirements  for  noncoherent  pro¬ 
cessing  are  two  or  three  orders  of  magnitude  less  than  for  coherent 
processing. 

•  Noncoherent  processing  is  much  more  robust  than  coherent  processing 
with  respect  to  relative  motion  of  the  target  and  the  radar. 

•  Coherent  processing  provides  target-to*clutter  contrast  that  is  superior 
to  noncoherent  processing. 

This  study  does  not  consider  the  role  of  wavelet  based  signal  analysis 
methods  for  noncoherent  radar  data  processing. 


D.2  Characteristics  of  the  Millimeter  Wave 
Radar  Data  used  in  this  Study 

Millimeter  wave  radar  data  was  gathered  on  a  test  range 
by  Martin  Marietta.  The  radar  return  signals  were  co¬ 
herently  detected  and  the  complex  data  recorded  into  3 
data  sets. 

1.  Calibration  radar  reflectors  which  consisted  of  an  array  of  two  dihedral 
reflectors  and  one  trihedral  reflector.  These  data  were  taken  with  every 
target  measurement,  but  only  50  records  were  provided  for  analysis 

2.  M60A3  Tank 
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Figure  D.l:  Real  Component  of  Trihedral  Calibration  Data 

3.  M35  Truck. 

Figures  D.l  to  D.6  provide  examples  of  the  raw  range  data  and  the  DFT  for 
the  calibration  reflectors.  Figures  D.7  to  D.9  are  examples  of  data  taken  of 
a  military  target  vehicle,  an  M60  tank  in  this  case. 

The  vehicles  were  located  atop  a  rotating  platform  adjusted  such  that  the 
depression  (look-down)  angle  was  15  degrees.  The  platform  rotated  at  about 
1  degree  per  second. 

The  data  sets  for  the  M35  truck  and  M60  tank  consist  of  14,000  records 
each.  A  record  consists  of  fotir  complex  valued  vectors  of  length  63.  Each 
vector  corresponds  to  a  single  set  of  63  transmitted  sinusoidal  waveforms 
(either  L  or  R  polarized).  The  first  and  second  vectors  in  each  record  cor¬ 
respond  to  a  L  polarized  transmission  followed  by  demodulation  with  a  left 
and  right  (respectively)  polarized  waveform.  The  third  and  fourth  vectors  in 
each  record  are  the  data  for  a  R  polarized  transmission.  Thus  each  record 
contains,  in  order,  polarization  modes  LL,  LR,  RL,  and  RR. 

The  time  between  consecutive  sets  of  transmitted  waveforms  (i.e.  a  pur 
of  consecutive  rows)  equab  6.3fi  seconds,  therefore,  each  data  set  represents 
176.4  degrees  of  rotation. 

The  frequencies  within  each  transmitted  set  span  a  bandwidth  of  £  =  500 
Megahertz  and  decrease  linearly  .  Thus  the  transmitted  waveform  for  the 
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Figure  D.4:  Real  Component  of  Dihedral  Calibration  Data 


Figure  D.5:  Imaginary  Component  of  Dihedral  Calibration  Data 
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k-tb  column  within  a  given  pair  of  rows  equals  sin(2T[/e  —  {k  —  62)£/63]0 
where  /c  equals  the  reference  frequency. 

For  both  vehicles,  the  angle  measured  from  head-on  (0  degrees)  coun¬ 
terclockwise,  as  a  function  of  record  index  k  is  given  by 

e{k)=^ak  +  b.  (D.12) 

For  the  tank  a  =  .0128  and  6  =  —23.9,  thus  the  angle  ranges  ^  from  —24 
degrees  to  156  degrees. 

For  the  truck  a  =  .0128  and  b  =  —34.6,  thus  the  angle  ranges  from  —35 
degrees  to  145  degrees. 


’The  data  was  divided  uto  *block8”  of  2000  teeords  for  ease  of  shipment  on  floppy 
disks.  It  was  discovered  late  in  the  study  that  one  "block”  of  tank  data  was  corrupt^. 
These  data  were  simply  skipped  leaving  12,000  records  of  M80  data. 
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Appendix  E 

Baseline  Target  Extent 
Calculation  Method 


This  section  details  the  baseline  method  for  extracting  target  extent  from  a 
phase  history  data  sequence  zt  < 


1.  compute  the  Discrete  Fourier  Transform  (DFT)  z*  of  the  complex  radar 
data, 

2.  compute  the  squared  modulus  of  the  transform  results 

3.  establish  a  threshold  T  >  0, 

4.  compute  an  estimate  of  the  noise  floor.  Since  the  data  provided  does  not 
have  any  ’‘empty’*  range  celb,  the  average  power  A  calculated  from  the 
16  range  cells  furthest  away  from  the  target,  was  used  as  the  estimate 
of  the  noise  floor.  This  was  the  best  estimate  available  since  the  true 
target  position  was  know  to  be  in  the  center  of  the  array  Zk.  The  first 
8  and  the  last  8  range  resolution  cells  were  used  to  estimate  A, 

5.  compute  the  indices  Ii  and  /a  of  the  first  and  last  entries  of  that 
exceed  T  x  A, 

6.  compute  target  extent  as  /]  —  /j  +  1. 
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